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ABSTRACT 

A careful study of the supersymmetric version of Pruisken's nonlinear 
a model for the integer quantum Hall effect is presented. The lattice 
regularized model is cast in Hamiltonian form by taking the anisotropic 
limit and interpreting the topological density as an alternating sum of 
Wess-Zumino terms. It is argued that the relevant large-scale physics of 
the model is preserved by projection of the quantum Hamiltonian on its 
sector of degenerate strong-coupling ground states. For values of the Hall 
conductivity close to e 2 /2h (mod e 2 /h), where a derealization transition 
occurs, this yields the Hamiltonian of a quantum superspin chain which 
is closely related to an anisotropic version of the Chalker-Coddington 
model. The relation implies that the ratio of magnetic length over po- 
tential correlation length is an irrelevant parameter at the transition. 
The superspin chain resembles a Id isotropic antiferromagnet with spin 
1/2. It has an alternating structure which however permits an invari- 
ance under translation by one site. The conductance coefficients of a 
quantum Hall system with iV small contacts translate into iV-superspin 
correlation functions which are governed by conformal invariance. The 
superspin formalism provides a framework for studying the crossover 
from classical to quantum percolation. It does not however encompass 
the frequency-dependent correlations of wave amplitudes at criticality. 
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1 Introduction 



When a two-dimensional (2d) electron gas is placed in a strong magnetic field B, the 
Hall conductance on is observed to be a rational multiple of e 2 /h near certain magic 
values of the density PJ. While Coulomb interactions are essential for the occurrence of 
fractional quantization, it is generally agreed that the observation of integer values an = 
ne 2 /h (n G Z) in low-mobility samples can be understood from a model of noninteracting 
electrons with potential disorder. It was first proposed in [|J that such a model undergoes a 
sequence of localization-delocalization transitions as the Fermi-energy Ep ~ 1/B is varied. 
Given the localization of all states in the Landau band tails, quantization of the Hall 
conductance has been explained by various theories, of which we mention the topological 
argument || and the scattering-theoretic approach [|J . What is less well understood is the 
critical phenomenon which occurs at energies E c close to a Landau level center where some 
electron states become delocalized. The picture that has emerged from experiments and 
numerical simulations for noninteracting systems ||[7j is that of some kind of second-order 
phase transition with a divergent correlation length £ ~ \E — E c \~ u , v = 2.3 ± 0.1. 

We begin with some general statements concerning conformal invariance at the integer 
quantum Hall derealization transition and other metal-insulator transitions in disordered 
electron systems, in Sect. 2. Two ways of regularizing the electron Green's functions are 
distinguished. The system may either be open with absorbing (i.e. conducting) boundary 
conditions at the leads or, alternatively, it can be taken to be closed but perturbed by an 
external field with finite frequency. Both an external field and absorption at a boundary 
cut off the long-time correlations and in this sense act as regulators. The introduction of 
a frequency u sets a length scale, L w , which breaks conformal invariance. Moreover, the 
amplitude of critical Green's functions diverges for L w — > oo (or uj —>■ 0). An interpretation 
of these Green's functions in the framework of conformal field theory || must involve pri- 
mary fields with scaling dimensions that are negative. The situation is somewhat different 
when the first option of regularizing by absorption at a boundary is utilized. We suggest 
the novel procedure of keeping the size and relative distance of the absorbing contacts 
finite while sending the system size to infinity. In this limit the d.c. conductance, or 
current-current correlation function, decays algebraically with increasing distance between 
the contacts, and we expect it to expand in conformal fields that have positive scaling 
dimensions only. 

To promote the reasoning of Sect. 2 to a quantitative theory with the possibility of 
getting exact analytical solutions, we must turn to field-theoretic techniques. An effective 
field theory for the integer quantum Hall effect was first formulated by Pruisken and 



collaborators PJICj- The field theory is of the nonlinear a model type and has the action 



S[Q] = \f dxdy (a° xx tv(dQ) 2 - a° xy tvQ[d x Q, d y Q\) . (1) 

The matrix-valued fields Q map the configuration space M. of the 2d electron gas into a 
coset space U(2n)/U(n) x U(n) with n = due to the use of the replica trick. By the 
derivation of the model, the bare coupling constants a xx and a xy are identified as the 
longitudinal and Hall conductivities, respectively, calculated in the self-consistent Born 
approximation and measured in units of e 2 /h. The field theory with action ([!]) attributes 
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to a xy the meaning of a topological angle 9 = 2na xy multiplying the topological density 
trQ[d x Q, d y Q]/8. When M. has no boundary, the topological density integrates to 2ixi 
times an integer and the physics of the model ([1]) is periodic in a xy with period 1. 

Pruisken's work on the model (|IJ) focussed mainly on the semiclassical instanton 
approximation [TJ,|12|, from which he proposed a two-parameter renormalization group 
flow with one (or possibly two) unstable fixed point(s) at a xy = 1/2 (mod 1) and some 
a xx = a xx- This (or one of these) fixed point (s) is conjectured to be the integer quantum 
Hall derealization fixed point seen in numerical and real experiments. While Pruisken's 
conjecture has great intuitive appeal and is consistent with experiments, it has not led 
to quantitative results, beyond predicting the critical point to be isolated - with no con- 
tinuous symmetry spontaneously broken and Goldstone bosons absent. In particular, the 
observed universality o* xx = 0.5 ±0.1 [|13| was not foreseen, and a computation of critical 
indices is still lacking. This is not surprising since the critical region is beyond the reach 
of the instanton approximation, which is valid for a xx ^> 1 while the integer quantum 
Hall (IQHE) transition occurs at a critical value a* x of order unity. Another, more serious 
difficulty arises from the constraints imposed by conformal invariance. The field theory 
(D at n = and a xy = 1/2 (mod 1) must be critical with a divergent correlation length 
- or else it would not describe the IQHE transition. We therefore expect it to flow under 
renormalization to a conformally invariant field theory H. Now, as was pointed out by 



Affleck [[nj, conformal invariance of a 2d field theory with continuous symmetry implies 
that there exist two sets of conserved currents, one for the right-moving and another one for 
the left-moving charges. Unfortunately, the nonlinear a-model ([]]) has just one manifestly 
conserved current, which is its Noether current. (Presumably, in order for the required 
separation into left and right movers to become visible, one has to restrict the a model to 
its low-energy limit.) What happens is well understood for the case of one replica (n = 1). 
Affleck and Haldane |Tj| have argued convincingly that the U(2)/U(l) x U(l) nonlinear 
a model (also known as the 0(3) model) at 9 = n is not scale-invariant but renormalizes 
to the SU(2) k = 1 WZW model - a conformal field theory with central charge c = 1 
and two independent conserved currents J{z) and J(z). This conclusion was confirmed by 
Shankar and Read ||16|| . By analogy, we expect Pruisken's model at 9 = n and n = to 
flow to some conformal field theory of the Wess-Zumino type. We are thus led to reject 
Pruisken's assertions about the fixed-point structure of the model (fjj at n = 0. This model 
has neither one nor two nontrivial fixed points at 9 = n, but has no such fixed points at all! 
(Note that we are not calling into question the two-parameter scaling theory ||17|| , which is 
arguably correct when proper meaning is given to the parameters a xx and a xy . What we 
are saying is that these parameters have no interpretation as running coupling constants 
of the model (|).) 

The lesson learned is that, if our goal is to approach the IQHE fixed point by analytic 
means, we had better abandon the critical (at 9 = rr) but scale-dependent field theory 
(|1]). This logical step was taken by Affleck jn| who exploited the relation of Pruisken's 
U(2n)/U(n) x U(n) replica field theory to certain SU(2n)-invariant quantum "spin" chains. 
Such a relation first appeared in the work of Haldane [[19| who suggested that antiferro- 
magnetic spin chains with SU(2)-invariance and large spin S map on the U(2)/U(l) x U(l) 
nonlinear a model with topological angle 9 = 2nS. Applying an inverse of this mapping 
jl^l to Pruisken's model for n = 1, one gets an interpretation of the Hall conductivity 
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a® y = 1/2 (mod 1) as the spin S = 1/2 of a Id isotropic antiferromagnet. Even more in- 
triguingly, the hypothesis of a sequence of localization-delocalization transitions occurring 
at a® y = ±1/2, ±3/2, ±5/2 etc. translates into the accepted statement that a Id isotropic 
S = 1/2 antiferromagnet, while being in a massive dimerized phase for a staggered site-site 
interaction, becomes critical as the staggering is turned off. 

Affleck's approach, though well conceived, ultimately failed. An early prediction [y = 
1/2, gotten from a further mapping of SU(2n) spin chains onto SU(2n) WZW models) was 
later withdrawn on the grounds that the phase transition at 9 = n is likely to be of first 
order for n > 2, making it difficult if not impossible to extrapolate to the replica limit 
n = 0. This did not come as a surprise. Problems with the replica a models of localization 
theory occur not only for the integer quantum Hall transition, but are notorious and show 
up whenever one tries to make calculations in a nonperturbative regime where observables 
are sensitive to the global topology of the field space. (The first example of a discrepancy 
with exact results was reported in ||20|| .) Here then is a summary of how I view the 
situation: i) The nonlinear a model ([!]) does not have the right symmetry structure to 
accommodate conformal invariance and is therefore a bad starting point for extracting the 
critical properties, ii) Replica field theory has failed. The replica limit is not smooth and 
one does not know how to extrapolate to n = 0. 

Clearly, a new analytical approach to the IQHE transition is very much called for. 
Such an approach was recently proposed by Ludwig et al. [^TJ who introduced a class of 
model systems which undergo a phase transition with a jump in the Hall conductivity even 
in the absence of disorder. This transition is in the Ising universality class and has v = 1. 
Adding various types of disorder to the pure system, Ludwig et al. found a rich array of 
fixed points and lines, but they were unable to control the perturbations that drive the 
system to the generic IQHE fixed point of interest. It remains to be seen whether any 
quantitative results for the proper IQHE transition will be obtained from this approach. 
In the present paper we will follow a more conservative line of attack. It is well known that, 
in addition to the replica trick, there exists another (and actually superior) method for 
calculating disorder averages of products of Green's functions for noninteracting electron 



systems. This is the supersymmetry technique pioneered by Efetov |22;1 , and it was used 
in [23] for the derivation of a supersymmetric analog of Pruisken's model (|l|). This latter 



model, which enjoys a mathematically sound foundation, is the object we shall study in 
the present paper. 

The supersymmetric model was rejected by Affleck on the grounds that its fermions 
do not satisfy the spin-statistics theorem. While this is indeed so, I fail to see why it should 
disqualify the model. To be sure, the supersymmetric model was not constructed so as 
to satisfy the narrow requirements of canonical quantum field theory! Rather, it arose 
from a careful study of noninteracting disordered electron systems, mapping the problem 
of calculating ensemble averages of products of single-electron Green's functions on the 
statistical mechanics problem of computing the Green's functions of a supersymmetric 
field theory. When viewed as a statistical mechanics model, the supersymmetric model 
makes perfect sense. In fact, it possesses a transfer matrix that passes for a well-defined 
object even by the most rigorous of mathematical standards. It is therefore natural to 
go back to Affleck's perceptive ideas [[18| and try and get them straightened out by using 
supersymmetry. This is what we shall attempt to do in this paper. Of course, things will 
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be rough going at times because we'll have to make do without the luxuries of canonical 
quantum field theory (Hilbert space, normalizable states, canonical quantization, unitary 
representations etc.) Nevertheless, as we shall see, the first half of Affleck's programme 

- mapping of the nonlinear a model on a quantum spin chain - can be executed. The 
second step - mapping of the spin chain on a field theory of the Wess-Zumino type and 
computation of the critical indices - is currently under investigation. 

In Sect. 3 we carefully pass from the 2d Euclidean supersymmetric field theory to its 
Id quantum Hamiltonian. Our treatment was inspired by a pedagogical paper of Shankar 
and Read 1TB] on the 0(3) nonlinear a model. Setting up the transfer matrix on a mesh of 
lattice points and taking the anisotropic limit, we obtain a Id lattice Hamiltonian which is 
a sum of two terms: a "kinetic energy" governing the field dynamics on individual lattice 
sites and a "potential energy" which couples fields on neighboring sites and enforces local 
order. The kinetic energy operator for a site located on one of the boundaries of the 
quantum Hall system is a supersymmetric generalization of the Hamiltonian of a charged 
particle moving in the field of a (fictitious) magnetic monopole, with the monopole charge 
being the integer m closest to 2cr° . We call this operator the "monopole Hamiltonian" 
for short. The value of the monopole charge is quantized by the condition that the Dirac 
string singularities of a corresponding magnetic vector potential be unobservable. For half- 
integral values of <r° , the single-site dynamics in the bulk of the system is governed by the 
monopole Hamiltonian with charge m = ±1, and the potential energy, which is staggered 
in general, becomes translationally invariant. 

The ground state of the monopole Hamiltonian is constructed in Sect. 4. The outcome 
of this calculation can be understood from the simple analogy to the Landau problem of 
a single 2d electron in a magnetic field. In a suitable gauge, the lowest Landau level is 
spanned by holomorphic functions of the complex coordinate z = x + iy, multiplied by a 
Gaussian. Reversal of the magnetic field turns holomorphic functions into antiholomorphic 
ones. Similarly, the lowest level of the monopole Hamiltonian of the a model is a space 
V (V*) of holomorphic (antiholomorphic) functions for monopole charge m > (to < 0). 
The Gaussian gets replaced by something else since the field space is not flat but curved. 

The significance of the integer to becomes particularly clear in Sect. 4.6. In previous 
work on Pruisken's nonlinear a model [ |i~T| , only the local (longitudinal and Hall) con- 
ductivities were calculated. However, it is clear that these are insufficient for a complete 
characterization of transport in a quantum Hall system - as in any mesoscopic conductor 

- at T = 0. What is needed is the entire set of conductance coefficients p3f . Computing 
these for a circular geometry with iV probes, we obtain a result which agrees with what is 
expected from phenomenological theory [Q: the nonvanishing conductance coefficients for 
even m are given by \m\ - the number of "edge states". 

The sign of the monopole charge m = ±1 alternates on sites in the bulk, resulting in an 
alternating sequence of single-site ground state spaces V ', V* , V, V* etc. Our calculation 
of Sect. 4.6 suggests that a state in V ® V * <E> V <8> V* <E> . .. can be thought of as representing 
a snapshot of the evolution in "time" (i.e. along a quasi-lcf strip) of the distribution 
of squared transfer matrix elements in an anisotropic version of the Chalker-Coddington 
model [0. Guided by this interpretation, we argue that the relevant large-scale physics 
of the a model is captured by simply projecting the site-site interaction on the space 
of strong-coupling ground states V <S> V* <g> V <S> V* <E> ... . The result of doing such a 
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projection is a quantum (super-)spin chain Hamiltonian specified in detail in Sect. 5.1. 
The possibilities for obtaining exact solutions of this or related Hamiltonians are under 
current investigation. 

We conclude this introduction with some comments on the style of presentation. This 
paper is intended to be fairly self-contained and complete, which explains its length. In 
Sects. 3-4, familiarity with super symmetry is needed only for some calculational details, 
many of which have been relegated to several appendices. The methods used in this 
paper are robust and apply to any symmetric space G/K, whether it is super or not. 
Indeed, everything we do is easily transcribed to the replica a model by simply replacing 
supertrace by ordinary trace, superdeterminant by ordinary determinant etc. and omitting 
the sign factors from superparity. We take advantage of this fact by illustrating many of 
the concepts at the familiar example of a two-sphere S 2 ~ U(2)/U(l) x U(l). 

The material of Sect. 2 is largely independent of the rest of the paper and will not be 
used until Sect. 5.1. 

2 Conformal Invariance of Critical Disordered Electron Sys- 
tems 

Conformal invariance, which occurs in models of field theory and statistical mechanics at 
a critical point [Q, is an especially powerful symmetry in two space dimensions, the group 
of conformal transformations being in this case infinite-dimensional. Conformal invariance 
connects correlation functions in different geometries with each other; for example, the 
conformal mapping z h> lnz relates the correlation functions in the plane to those on a 
cylinder. Relations of this kind are being widely used for the computation, numerical and 
otherwise, of critical indices. 

This second section of the paper is concerned with the question whether conformal 
field theory has applications to critical wave propagation in 2d disordered media. It is by 
now well known |25| that the interference of waves multiply scattered by strong disorder 



leads to localization, if the coherence volume is large enough. The points in parameter 
space where the transition from diffusive motion to localization occurs, are characterized 
by a diverging correlation length. It is therefore reasonable to ask, and the question has 
been asked before [[26] , whether conformal field theory makes statements about the critical 
correlations at such a point. To be specific, let us consider in the sequel the motion of a 
single electron in a disordered solid, although the reasoning used may have generalizations 
to other systems. Then there are (at least) two universality classes where a critical point 
occurs in two dimensions: (i) systems in a strong magnetic field (quantum Hall systems), 
and (ii) time-reversal invariant systems with spin-orbit scattering. 

For a system of class (i) or (ii), we take (■) to mean the average over some statistical 
ensemble of Hamiltonians H which is invariant under the group of Euclidean motions of 
the infinite plane. Let 5(E) = n~ 1 e/(e 2 + E 2 ) be a regularization of the 5-distribution, 
and let S E (x,y) denote the matrix element of 5(E — H) between the eigenstates of the 
position operator at x and y. With these definitions we consider the averages 

(S E (x,x) p S E +hLo(y,y) q ) (p,geN) and (S E (x,y)S E+hul (y,x)). (2) 

The first of these measures the correlations of the local density of states, while the Fourier 
transform of the second one yields the dynamic structure factor, which can be parametrized 
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in terms of a wave- vector and frequency-dependent diffusion constant P^j . By the spectral 



expansion Se(x, y) = ^2 ^i{x)5{E—Ei)ipi{y)^ the correlation functions (||) can expressed in 
terms of the eigenfunctions of H with energies E and E + hui. Introducing the symbol ((•)) 
to to denote ensemble averaging in combination with energy averaging by the regularized 
5-distribution, we can write the correlators (||) in the short-hand form 

K AB (x-y ] E,u) = ((A(x)B(y))) (3) 

where A(x) (B(y)) is some polynomial in the eigenfunctions and their complex conjugates 
evaluated at the point x (y). Kab is translationally and rotationally invariant, by the 
Euclidean invariance of the probability measure defining (•). 

At a critical point, Kab(x — y) is expected to show algebraic decay ~ \x — y\~ v 
with some positive exponent rj. When the infinite plane is replaced by a cylinder with 
circumference L, finite-size scaling theory states that Kab should decay exponentially 
~ exp(— \x — y\/XL). One might now think that, by the usual hypothesis of conformal 
invariance at a critical point, 77 and A are related by 

77 = 1/Xtt. (4) 

This however is false and we will review the reason why in Sect. 2.1, by giving a simplified 



presentation of the work of Wegner |f28[ . The lesson learnt from this will prompt us to 
point out in Sect. 2.2 a different class of correlation functions where the relation (|]) does 
hold. 

2.1 Finite frequency breaks conformal invariance 

The message of this subsection is easy to state. In addition to the correlation length £ 
there exists a second length scale, L u , which is set by the frequency ui. The correlation 
functions Kab, Eq. ©, re tain their dependence on at the critical point, even in the 
limit — > 00. The precise form of these correlation functions has been calculated by 
Wegner [^8| in some cases. We will now repeat Wegner's basic argument. 

Let be defined as the linear size of a system with level spacing frw: L w = {fijjjp)~ 1 / 2 , 
where p is the density of states per unit of energy and area. At the critical energy E c where 
the correlation length £ ~ \E — E c \~ v diverges, Kab takes the form 

K A b(z, E c ,u) = k AB (\z\,Lu) (z = x-y). (5) 

A necessary condition in order for arguments from conformal invariance to apply is that 
k,AB should be ignorant of the length scale L u . For L w large and fixed, let us therefore 
specialize to the limit \z\ <C L u where we expect power law behavior 

k A B{\z\,L u ) ~co\z\-^ 

with some exponent t\ab- 

We can obtain more information about t\ab and the dependence of cq on by using 
the real-space renormalization group (RG) transformation for disordered single-particle 
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systems introduced by Wegner |2l|. To keep the technicalities to a minimum, we assume 
that A and B are scaling fields {i.e. eigenfields of the RG transformation) with scaling 
dimensions and Ag, respectively. A change of cutoff scale a i— > b±a then results in 

A(x) i- b^Aix/h), B(y) i ► b^ AB B(y/h). 

Choosing b\ = \x — y\/a = \z\/a, we get 

k AB (\z\,L w ) = (\z\/a)^ AA+A ^k AB (a, L u a/\z\). 

Note the dependence of the right-hand side on \z\ through both the algebraic prefactor 
and the second argument of k A B, which is L u a/\z\. To extract the full | z (-dependence 
of the correlation function (|^), we must therefore perform a second RG transformation 
bid i — ► 62^1^, taking L u a/\z\ to some fixed value L^a/ (\z\b2) = L . In calculating the 
effect of this transformation, we may no longer renormalize A and B separately, since the 
RG transformation a ^ b±a = \z\ has already brought the points at which A and B are 
evaluated within a distance of one unit of the cutoff scale. Instead, we must renormalize 
the local object A(x /b\)B(y /b\) ~ (AB)(x/bi). The product AB will not in general be a 
scaling field, but it can be expanded in terms of such fields. If A A b denotes the scaling 
dimension of the most relevant field occurring in this expansion, (AB)(x) is renormalized 
under a 1— > 62 a by 

(AB)(x) 1 — > &2~ Aab (AB)(x/b2) + subdominant terms. 
Taking 62 = L u a/(Lo\z\) = const x L w /\z\, we finally obtain: 

k AB (\z\,L u ) ~ | Z |-A.-A B +A ABL -A AB? (6) 



which expresses the dependence of k A B on both \z\ and by three scaling dimensions 
A A , A^, and A^^- On physical grounds, the correlation function k A B must neither 
increase with distance \z\ nor decrease with L u , uj ~ L~ 2 being a distance in energy 
space. Therefore, A A b < and A A b < A A + Ab- For the case of potential scattering 
(orthogonal universality class), Wegner has shown the scaling dimensions A A and Ab for 
the first choice in (|]) to be zero or negative in 2 + e dimensions. (Moreover, the scaling 
dimensions decrease without bound as the powers p and q increase.) We expect this to be 
true in general. Thus, the scaling dimensions Ao (O = A, B, or AB) are subject to the 
inequalities 

A o <0, -A AB >-A A -A B . 

The second of these can be understood independently of @ by observing that the limit 
L w — > 00 restricts the contributions to k A s to those from a single eigenfunction i/j with 
energy close to E c . The inequality is a manifestation of the existence of large amplitude 



fluctuations of such an eigenfunction, indicating its multifractal nature [30 |. 

We now return to the question concerning the relation of t\ A b = A A + Ag - A^# to 
the Lyapunov exponent for exponential decay of k A s on the cylinder. The derivation of 
relation (|]) assumes that A and B are (spinless) primary fields which are dual to each other 
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and belong to a unitary representation of the Virasoro algebra, so that A a = = A > 
and (A(x)B(y)) ~ \x — y\~ 2A . Clearly, these assumptions are violated in the present 
context. We are therefore ill-advised [|26| to expect a relation of the form (^) to hold. (Of 



course, this does not yet rule out the possibility of some generalization of ([|) still being 
valid.) 

Could there exist some sort of sensible, nonunitary conformal field theory that has 
an infinite number of primary fields with negative scaling dimensions to accommodate the 
correlation functions (^|)? Whatever the answer to this question may turn out to be, it 
will be constrained by the observation |28[] that ((A(x))) for A(x) = Se(x,x) p and p>2 
diverges as the regularizing imaginary frequency is — > 0. This divergence is suggestive of a 
field theory with unbroken noncompact symmetry - just as (A) = for an order parameter 
field A would signal unbroken compact symmetry. In fact, ever since the seminal work of 



Schafer and Wegner [31] and of Efetov []22] it has been understood that the problem of 
calculating disorder averages of powers of wave amplitudes maps on field theories of the 
nonlinear a model type with noncompact symmetry. We will study such a field theory in 
Sects. 3-4. The nonpositivity of the dimensions of a subset of its scaling fields turns out to 
be intimately related to the fact that these fields are improper, i.e. their application to the 
vacuum produces states that are non-normalizable in the absence of external frequency as 
a regulator for non-compactness. 

2.2 Conformal invariance and d.c. conductance 

We consider again a 2d noninteracting disordered electron gas, either in a strong magnetic 
field or with spin-orbit scattering. The configuration space M. for one electron may be 
finite or infinite. An important difference to earlier is that we now envisage the iV-terminal 
geometry shown in Figure 1. There are iV contact areas where electrons are exchanged 
between M. and the leads. Applying a small voltage (or, more precisely, a small shift of 
chemical potential) to lead q we get a small current flowing in lead p. The constant of 
proportionality is called the d.c. conductance, or conductance coefficient, g pq . From linear 
response theory one can derive a formula for g pq as a current- current correlation function 
2^| : g pq ~ Tv(v p G + v q G~) in symbolic notation where v p is the component of the velocity 
operator normal to the contact area C p , and the Green's functions G r± are those of the 
total system (including the leads) evaluated at the Fermi energy Ep. It turns out that an 
alternative formulation on reduced configuration space Ai (without the leads) is possible. 
If r is a certain nonlocal kernel determined by the lead geometry, it has been shown [j32[ 
that 

9 pq = / dw I dx dy dz T(w : x)G + (x : y)T(y : z)G~(z,w). (7) 

J Cp J C p J C ' q J Cq 

G ± are now the Green's functions of the isolated system subject to certain boundary 
conditions on the contact areas. Qualitatively, G^ = (Ep — H\ M ± zr) _1 . The exact 
formulation has been given in fl3"2fl . (Note that Eq. (8) of |3^] is misprinted: G~ should 



read G~^ instead.) Formula ([F[) is exact for noninteracting electrons at zero temperature 
and for a magnetic field of arbitrary strength. We are measuring conductance in its natural 
units e 2 /h. 

The energy argument of the reduced Green's functions G^ in (|7]) is a real number 
E = Ep. Why are we permitted to set to zero the infinitesimal imaginary part ±ie which 
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is usually added to energy to make Green's functions well-defined? The answer is that the 
boundary conditions satisfied by the reduced Green's functions spell out the "openness" 
of the system, i.e. the possibility for electrons to leave M. at a contact area C r and enter 
the lead r. The damping caused by the escape of probability from Ai shifts the spectrum 
of H | M off the real energy axis and makes the limit e — > exist. This is the crucial point: 
conductance - as opposed to the Green's functions for an isolated system - does not require 
any regularizing real or imaginary frequency in order to be well-defined. 

Let us now take the disorder average (g pq ) and set energy right to its critical value 
E c where the correlation length diverges. When the spatial extension of the contact areas 
is comparable to the system size, as is the case in most experimental geometries, the 
average conductance will become some scale-invariant number (g pq ) = g* q . We find it 
more interesting however to consider the opposite case of small contacts, see Figure 1. A 
procedure that makes perfect sense is to take the thermodynamic limit of infinite system 
size while keeping the relative distance and size of the contacts fixed. We expect the 
conductance coefficients to be well-behaved and finite in this limit. With system size 
having gone to infinity, there do not remain any scales over which to measure correlations. 
(Of course, one can still assign some length scale to the "strength" or "decay width" |33 



of each contact. However, since the contacts are no more than a local perturbation on 
the system, they do not set any scale in the sense of scaling theory. Note that, in spite of 
this, they remain effective as regulators for conductance even in the thermodynamic limit.) 
With all scales absent, the average conductance (g pq ), and all moments of conductance, 
must decay algebraically with increasing distance between the contacts. (This will remain 
true even for large contacts, provided that their size is kept fixed as the distance between 
them is sent to infinity.) 

Within the framework of the effective field theory described in Sect. 3, we can rephrase 
our assertion as follows. If the linear extension of the contacts is less than or comparable to 
the microscopic cutoff scale of the field theory, which is set either by the magnetic length or 
by the elastic mean free path, the contacts can be regarded as point-like. The conductance 
coefficients (g pq ) for iV terminals then translate into iV-point functions 

(O p (x p )O q (x q )Y[O r (x r )), (8) 

r 

with local operators Oi(x), of the field theory. On physical grounds, these iV-point functions 
must decay algebraically on all but the shortest distances when E = E c . There aren't any 
length scales left and, moreover, the disease of divergent amplitudes is now absent. It is 
therefore reasonable to expect that the iV-point functions (JD can be expanded in terms 
of scaling fields with positive dimensions of an underlying conformal field theory. (As a 
matter of fact, we shall see in Sect. 5.1 that the iV-point functions (|8|) restrict to dynamical 
iV-spin correlation functions of a superspin chain Hamiltonian for quantum Hall systems.) 
More precisely, because of the noncompact symmetry of the field theory, we expect a 
continuum of primary fields - with an associated continuum of scaling dimensions - 



to occur in the expansion of Oi(x). 

In summary, what we are saying is the following. Critical disordered electron systems 
possess a number of interesting correlation functions whose definition involves two (or 
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more) energies E + and E_. The effect of frequency uj = (E + — E_)/h is similar to the 
effect of an external magnetic field on a critical spin system: it sets a length scale and breaks 
conformal invariance. While conformal invariance is restored for the spin system when the 
external field is turned off, the same is not guaranteed to happen for the disordered electron 
system. Many of the correlation functions of interest diverge as uj — > 0. For example, the 
dynamic structure factor does. (This is a consequence of critical long-time correlations 
at the metal-insulator transition.) Any consistent conformal field theory interpretation 
of these correlation functions will have to confront the task of defining operator product 
expansion for an infinite set of primary fields whose spectrum of negative dimensions is 
unbounded from below. The situation is much better for correlation functions that can 
be defined at a single energy E = E c . An example is the d.c. conductance - or current- 
current correlation function - of an open system. While conformal invariance will govern 
the geometric variation of conductance for any contact size, things become particularly 
simple for point contacts: in this case the conductance coefficients for an iV-terminal 
geometry map on iV-point functions for local fields of an effective field theory. These 
iV-point functions are expected to have a conformal field theory representation involving 
primary fields with positive scaling dimensions only. 

I have included the material of Sects. 2.1-2 in this paper to motivate why the quest 
for the conformal limit of the 2d nonlinear a model holds more promise than just an 
explanation of the critical indices v = 2.3±0.1 Vi Tr = — 0.40±0.04 |J5j for the quantum 



Hall effect and v = 2.75 ±0.15 for systems with spin-orbit scattering [^6|. Conformal field 
theory will also predict interesting algebraic decay of the conductance and its moments for 
an iV-terminal geometry with small contacts. 

3 Hamiltonian limit of the nonlinear a model 

In the remainder of this paper we will make a fresh field-theoretic attack on the integer 



quantum Hall effect, starting from a supersymmetric version [23[ of Pruisken's nonlinear 
a model |J. Most of what we do in Sects. 3-5 is general and can be transcribed to systems 
with spin-orbit scattering without difficulty. I hope to return to such systems in a future 
publication. 

We start out by introducing the nonlinear a model from a gauge theory point of 
view in Sects. 3.1-2. We are not aiming at a complete solution of this model here. All 
we care about is its low-energy limit, describing the long wave length physics of integer 
quantum Hall systems. This is important to keep in mind, since restriction to the low- 
energy domain gives us much freedom to manipulate the field theory. We use this freedom 
in Sects. 3.3-5, where a careful derivation of the quantum Hamiltonian of the a model is 
presented. (The well-known connection between Feynman's path integral and Schrodinger 
wave mechanics allows us to relate the 2d classical field theory to a Id quantum field 
theory.) The possibility of such a Hamiltonian formulation was first discussed briefly in 



[10]. Our approach follows the work of Shankar and Read |16|] on the 0(3) nonlinear a 
model with topological angle 6 = n to some extent. The essential idea is to rewrite the 
topological density as an alternating sum of Wess-Zumino terms. Gauge ambiguity of the 
Wess-Zumino term quantizes the values of the topological coupling constant for which such 
a rewriting can be done. The 0(3) model is somewhat special in that its field space is a 
two-sphere and the topological density has an intuitive meaning as (the pullback of) the 
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solid-angle two-form. The simplifications resulting from this allow Shankar and Read to 
take some short cuts that are not available to us here. 



3.1 Review of Pruisken's model (supersymmetric version) 

An effective field theory for noninteracting electrons moving in two dimensions under the 
influence of a strong magnetic field and a Gaussian white noise potential, was first derived 
by Pruisken ||. He pointed out that a strong magnetic field has the effect of adding a 
so-called topological term to the Lagrangian of the nonlinear a model for noninteracting 



disordered electron systems developed by Wegner, Efetov, and others [|37| , |31| , P2| . A closely 
related model with internal supersymmetry, avoiding the replica trick used in PJIIJ, was 



subsequently derived by Weidenmiiller [g3| . In the present paper we shall study the latter 
model, extended by the addition of an extra pair of indices to the a model matrix field. 
Such an extension is necessary for the calculation of disorder averages of products of n 
retarded and n advanced Green's functions (n > 1), which is what is needed when one 
wants to go beyond the average of the conductance and calculate its variance and higher 



cumulants as well ||38|| . (Another benefit from making the extension is that it forces us 
to use the generic structures of the theory.) We will not give yet another derivation of 
the a model here but refer to the original papers for this |3l] , ^2| ,|9| j2^1 . The nonlinear a 
model is usually presented in terms of matrix fields Q satisfying the nonlinear constraint 
Q 2 — 1. This presentation, while adequate for some purposes, is rather useless for most of 
what we intend to do in this paper. So, rather than wasting time by first writing down the 
Q-field formulation and then changing to something better, we shall present the "good" 
formulation and the proper structures of the model right away. 

We start from the associative algebra of An x An supermatrices acting on the tensor 
product of three linear spaces: (i) superspace or Boson-Fermion (BF) space, (ii) Advanced- 
Retarded (AR) space and (iii) replica or Extra (E) space. The first two have dimension 2 
and the last one dimension n. Invertible An x An supermatrices g satisfying the condition 

9~ 1 = V9 t V with n = ((E n ) BF <8> (a 3 ) AR + (E 22 )bf ® Iar) ® (9) 

form a noncompact Lie supergroup G = U(n, n\2n). (Here as always {i2ij}i ) j=i,...,jv is the 
canonical basis of End(C JV ) and o\ = E\2 + E21, 02 = —%E\2 + 1E21, 03 = E\\ — E22 for 
N = 2 are the Pauli matrices.) The subgroup of elements in G that commute with 

A = l B F ® (0-3) AR ® 1e 

is denoted by K = U(n|n) x U(n|n). G is projected onto G/K, the space of left cosets 
gK (g G G), by g 1— >■ 7r(g) := gK. It is this coset space G/K which is the field space 
of the nonlinear a model. The base manifold (or "body" in the language of deWitt ||39|1 ) 
of G/K is the direct product of U(n,n)/U(n) x U(n) with U(2n)/U(n) x U(n). Both 
of these are Riemannian symmetric spaces [0, the former of the noncompact and the 
latter of the compact type. We sometimes refer to them as BB (Boson-Boson) and FF 
(Fermion-Fermion) space, respectively. 

To write down the a model action - and at the same time prepare later analysis - we 
identify various structures on Q = Lie(G), the Lie algebra of G. Defining orthogonality by 
the natural symmetric form 

B(X,Y) = STrAy/2 , (10) 
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where STr is the supertrace, we denote by V the orthogonal complement of /C = Lie(K) 
in Q. The orthogonal decomposition of an element X E G hy Q = K, + V is written 
X = X K + X v . (The explicit formulas are X K = (X + AX A) /2 and X v = (X - AX A) /2.) 
With these definitions the principal part of the a model Lagrangian can be written in the 
form 

L (g, dg) = B ((g-%g) T , (g'^gh) . (11) 

(Repeated Greek indices are implicitly summed over the values \i = 0,1.) Pruisken's 
Lagrangian |§ contains in addition to ([Tl]) another term, the so-called topological density 
Ltop, which is formulated as follows. Let a linear operator J : V — > V be defined by 
JX = i[A, X]/2. J squares to minus the identity and preserves B. With its help we define 
Q(X, Y) for X, Y £ V by Q(X, Y) = B(JX, Y). This two-linear form is skew symmetric: 

n(Y,x) = b(x, jy) = b(jx, j 2 y) = -n(x, r), 

owing to the symmetry of B and the properties of J. It induces the topological density by 

W0> dg) = e^fi ((g-'d^r, {g^d^v) (12) 

where e^ v is the antisymmetric tensor. Both Lq and L top are invariant under a local change 
of representative g(x) i— > g(x)k(x) (k(x) £ K) and are therefore well-defined as Lagrangians 
on the coset field space G/K. In addition, they are invariant under global transformations 
g(x) I— > gog(x) (go £ G). Eqs. (JTj) and (|12D follow from the Q-field expressions L = 
- STr(<9Q) 2 /8 and L top = STr Qd^Qd„Q/8i on making the substitution Q = gAg' 1 . 

After these preparations, we introduce the nonlinear a model for an integer quantum 
Hall system with longitudinal conductivity a xx and Hall conductivity a xy by the functional 
integral 

Z = V[g] exp- / d 2 x (a xx L + ia xy L tov ) ■ (13) 
J J M 

o xx and a xy are measured in their natural units e 2 /h. M. is the configuration space 
of the two-dimensional electron gas, and the functional integration measure T>[g] equals 
Ox dgi<(x) where dgx is the uniform super-integration measure for G/K. (Note that the 
derivation of the nonlinear a model from a Gaussian white noise potential yields for a xx 
and a xy their bare (or SCBA) values, denoted by a xx and a xy in |||10|]. We here imagine 
that these coupling constants have already been renormalized and simplify the notation 
by dropping the superscript 0.) 

The functional integral ( |13D calls for some additional explanation. First of all, the 
invariant symmetric form B, Eq. (|T0|), has the crucial property of being positive definite 
on the ordinary part (or body) of V. This means that G/K is Riemann in its natural 
geometry induced by B. Second, an experimental quantum Hall system has a boundary 
dM., which is composed of insulating and conducting parts, giving rise to corresponding 
boundary conditions for the a model field in flT3|). What is particularly important is 
that the conducting boundary conditions cut off the spatially uniform fluctuations of the 
noncompact a model field - which would otherwise cause a divergence - by breaking the 
global G-symmetry of ( |13D down to a global K-symmetry. (The microscopic origin of the 
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damping is the absorptive term iT in the reduced Green's functions G ± = (Ep — H\ M ± 
zr) _1 .) This, in conjunction with the Riemannian nature of G/K, makes the functional 
integral (|13|) well-defined and convergent after UV regularization. Third, for the sequel we 
adopt the Corbino geometry M. = [Ri, R2] x S 1 , with two or more contacts placed along 
the boundary dM.. Alternatively, for a particularly neat theoretical model, we might 
imagine M. to be a two-torus S 1 x S 1 , with several point contacts attached to it to cut off 
the divergent a model fluctuations. What is important for our derivation of the quantum 
Hamiltonian in Sect. 3.5 is that one of the two directions on M. be periodic. This direction 
~ S 1 will be assigned the role of imaginary time. Finally, we remark that the functional 
integral Z, Eq. (|T3|), equals unity by supersymmetry [12] which serves to cancel all vacuum 
graphs. (This is true even in the presence of a boundary since the boundary conditions 
preserve K-supersymmetry.) To obtain the Green's functions of the nonlinear a model one 
adds sources to the Lagrangian and takes derivatives as usual. 

For a xy = 0, (|T3| ) is the nonlinear a model describing the long wave length excitations 
of noninteracting disordered 2d electron systems with time-reversal invariance broken by 
a weak magnetic field (unitary universality class). Perturbative renormalization group 
analysis p7|j22]1 shows its coupling constant l/cr xx to be driven to strong coupling at large 
scales. This has led to the widely accepted conjecture that all stationary single-particle 
states of such systems are localized in d = 2. At the same time, the plateau-to-plateau 
transition in integer quantum Hall systems cannot be understood if all bulk electronic 
states are localized. There is ample numerical and experimental || evidence that 
states do in fact become delocalized near the center of each Landau band. (The data points 
to derealization occurring only at a single energy.) From accurate numerical studies for 
the lowest Landau level, we know the correlation length exponent of this derealization 
phase transition to be v = 7/3 (or close to that value). Hence, if ( |i~3l) is to describe 
the underlying physics, the addition of the topological density L top must cause critical 
behavior, i.e. vanishing masses for some of the a model fields, at certain values of a xy . At 
first sight, though, L top appears to be an unlikely candidate for the job it is supposed to 
do. Ltop does not change the classical equations of motion, for it can be written locally as 
a total derivative. This is closely related to the fact that, for suitable boundary conditions, 
/ d 2 x L top is topologically quantized in integer multiples of 2ir [10,12], making it invisible 
in weak-coupling perturbation theory. Nevertheless, Pruisken and collaborators [10] have 
argued emphatically that the presence of the topological density does lead to decisive 
changes in the large-scale physics of the replica analog of the nonlinear a model (|T3"D. 
One of their arguments uses twisted boundary conditions and ideas stimulated by t'Hooft 
duality to show that some kind of phase transition must indeed occur somewhere in the 
range from a xy = m to a xy = m + 1. However, the precise nature of the phase transition 
is not specified by the duality argument. It could, in principle, be of first order with 



a finite correlation length. (This is in fact the conclusion reached by Affleck [|18| : the 
U(2n)/U(n) x U(n) nonlinear a model for n > 2 has a first-order transition at 9 = n.) 
As was mentioned in the introduction, the current algebra of the model ([13]) fails to meet 
the strict requirements set by conformal invariance in d = 2. By implication, the model 
( P"3"| ) cannot possess any nontrivial RG fixed points, contrary to what has been asserted by 
Pruisken et al. [10,11]. Rather, what must happen is that the model at a xy = 1/2 (mod 1) 
is driven under renormalization to the fixed point of some other field theory, which is not 
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known at present but is likely to be of the Wess-Zumino type. To summarize the situation, 
a sceptic will say that fifteen years after the discovery of the integer quantum Hall effect, 
one of its crucial features, the occurrence of a localization-delocalization transition in the 
center of a Landau band, is still lacking an analytic description. This very unsatisfactory 
state of affairs has been the motivation for the efforts reported in the present paper. 

3.2 Gauge fixing 

The nonlinear a model for noninteracting disordered electron systems is usually presented 
as a theory of fields Q taking values in some nonlinear space isomorphic to a coset space 
G/K. In Sect. 3.1 we gave its alternative presentation [10] as a theory of G-valued fields 
g with an invariance under local gauge transformations g(x) i— > g(x)k(x) for k(x) G K. In 
other words, we are here viewing the nonlinear a model as a kind of gauge theory with the 
gauge group K acting on the field space G by right multiplication. (However, as the gauge 
field has no dynamics, the nonlinear a model is not a gauge theory in the strict sense of 
the word.) The practical implementation of the gauge theory point of view requires gauge 
fixing, i.e. some prescription for separating the unphysical (or gauge) degrees of freedom 
from the physical ones. Gauge fixing in the present context is best described as follows. We 
choose some (locally) smooth map s : G/K — > G. The only requirement on s is that, with 
7r being the canonical projection G — > G/K, s has to satisfy tvos = identity. Equivalently, 
for every g G G the element k(g) defined by the equation 

g = s(n(g))k(g) (14) 

must be an element of K. Technically speaking, an s with this property is called a section 
of the bundle n : G — > G/K. For the purpose of illustration, consider the example 
G = SU(2) and K = U(l) (generated by ia^). The coset space SU(2)/U(1) is isomorphic 
to the two-sphere S 2 by the map #U(1) i— > go~%g~ l = Yl^=i n i°~i-> assigning to every coset 
fl'U(l) a unit vector n = (m, ri2, n^). If we parametrize S 2 by a polar angle 9 and an 
azimuthal angle <fi in the usual way, two examples of a section s : S 2 —>■ SU(2) are 

8 = e #<X3/2 e i0<7 2 /2 and s = e i4>a 3 /2 e iea 2 /2 e -i4>a 3 /2_ ^ 

The property n o s = id is easy to check using the isomorphism SU(2)/U(1) — S 2 . 

After this example we continue our general considerations. By definition, a section 
s assigns to a coset %(g) = gK one representative s(n(g)) G G. Since this assignment is 
(locally) one-to-one, s distinguishes some submanifold of G which is (locally) isomorphic 
to G/K. Thus, making a choice of s fixes the gauge: s identifies the subspace of physical 
degrees of freedom on the field space G, which are integrated over in the functional integral, 
while k : G — ► K defined by Eq. (|14[) accounts for the unphysical gauge degrees of freedom 
(not integrated over). There is of course much arbitrariness in the choice of s. However, 
given any two sections s and si we can always find a map k\ : G/K — > K such that s and 
Si are related by a gauge transformation s = siki. Thus, the arbitrariness of s precisely 
reflects the freedom we have in choosing the gauge. 

We were careful to insert the qualification "locally" in various places of the above 
discussion. As a matter of fact, s cannot be defined globally when G/K is topologically 
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nontrivial. To see this by way of example, let us focus again on SU(2)/U(1) ~ S 2 and 
consider the first of the two sections given in (|i~5|) . Since the azimuthal angle cf> is ill- 
defined on the north and south pole, s is singular at these points. (It is perhaps helpful to 
mention that a choice of s specifies, by formula (27) below, a choice of vector potential A 
for a magnetic monopole placed at the center of S 2 . The singularities of s then translate 
into Dirac string singularities of A.) We can move the singularities of s around at will by 
a gauge transformation, but a basic result of fibre bundle theory tells us that we can 



never make them disappear altogether. The same kind of topological obstruction occurs 
for our quantum Hall model space G/K, due to the presence of the topologically nontrivial 
space U(2n) /U(n) x U(n) in its base manifold. (In this case however the singularities are 
not points but are submanifolds of codimension two on U(2n)/U(n) x U(n).) 

We need to make a small addendum concerning the most general choice of section 
s : G/K — > G that will be permitted in this paper. It so turns out that one crucial step - 
going from Eq. (25) to Eq. (28) below - requires the values taken by s to be restricted to 
Go = SU(n, n\2n), the subgroup of elements of G with superdeterminant equal to unity. 
There are good reasons for not making the restriction to Go throughout our development. 
(Working with Go instead of G would cause a certain amount of inconvenience in Sect. 5.) 
Thus, we will base all our considerations on the pair G and K, except that a section s of 
the bundle G — ► G/K will always have to be a map s : G/K — ► Go- 

3.3 Quantum Hamiltonian for a xy = 

To make a continuum field theory such as ( |i~3l) well-defined, one needs to specify some 
prescription for regularization of the functional integral in the ultraviolet. Elementary 
particle physics has come up with quite a number of different regularization schemes for us 
to use. The approach we adopt here is to reformulate the two-dimensional Euclidean field 
theory as a one-dimensional quantum theory defined on a lattice in space. This particular 
regularization scheme has been used successfully by Shankar and Read Jl6j in their frontal 
assault on the 0(3) nonlinear a model with topological angle 9 = n. It is the optimal choice 
in the present context, too, as it will allow us to clearly identify the relevant low-energy 
degrees of freedom. 

The quantum Hamiltonian of a canonical 2d Euclidean field theory can be constructed 
in the following standard way. First, one re-interprets the Euclidean theory as a Lorentz- 
invariant field theory in 1+1 dimensional Minkowski space-time. Then, the Lagrangian of 
the latter is converted into a Hamiltonian function by carrying out a Legendre transform. 
And finally, the classical Hamiltonian is turned into a quantum operator by the process of 
canonical quantization. What is attractive about this procedure is that it avoids detailed 
manipulations of the functional integral. It turns out, however, that the supersymmetric 
field theory ( p~3|) does not have a Hilbert space for its space of quantum states, so it is not 
a quantum field theory of the canonical type. Canonical quantization, though applicable 
as a purely formal procedure, is therefore not really appropriate here. (The noncanonicity 
is probably one of the reasons why the model (|13D has not yet been embraced by the 
community at large but has been studied by a rather small group of people.) I emphasize 
that there is nothing ill-defined or suspect about our noncanonical field theory. After 
lattice regularization, we are faced with the bona fide statistical mechanics problem of 
computing a finite number of well-defined and convergent superintegrals. 
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The foregoing discussion motivates the no-nonsense approach we will adopt. Having 
UV-regularized the field theory by placing its fields on a rectangular lattice, we will take 
the anisotropic (or Hamiltonian) limit and use no more - and no less - than straightforward 
manipulations of the integrals over field variables to derive the quantum Hamiltonian. More 
precisely, what we will do is to carry out the familiar steps leading from Feynman's path 
integral (in imaginary time) to the Schrodinger picture of quantum mechanics. Aside from 
getting the proper structures into play, this approach has the added advantage of leading 
directly and without further ado to the expression for the quantum Hamiltonian optimally 
suited for later analysis. There exists one delicate aspect however: the topological term 
needs to be handled with care. To avoid dealing with too many issues at once, we here set 
a xy = and postpone the inclusion of the topological term to Sect. 3.5. 

As is usual in statistical mechanics, we define the quantum Hamiltonian H as (minus) 
the logarithm of the transfer matrix T connecting two neighboring rows of lattice sites, 
both at constant time. (The first of the two directions on M. = I x S 1 is decreed to be 
"space" I = [J?i, R2], the other "time" or, rather, imaginary time.) In order for H = — InT 
to come out to be a local and tractable operator, we need to make the system anisotropic by 
lowering the coupling, or "temperature" , in time direction while strengthening it in space 
direction. Such an ad hoc modification of the system is permissible if we are concerned (as 
we are) with the study of critical or low-energy properties only, since spatial anisotropy is 
an irrelevant perturbation at a critical point. Imagining the anisotropy to be caused by 
our use of two different lattice constants - a small one, ao, in time direction and a larger 
one, a±, in space direction - we make the replacement 

/ dxdr L (g : dg) — ► V" ( — A(g n , g n + eo ) + — A(# n , tfn+ej ) (16) 
Jixs 1 „ V a o ai J 

with the "distance" A(g,h) between g and h defined by A(g,h) = STr g~ 1 hAh~ 1 gA/4. 
The sum runs over all lattice points n, and eo and ei are the lattice vectors in time 
and space direction. Periodic boundary conditions g n +N e = 9n are assumed. Our dis- 
cretization scheme has the right naive continuum limit, as is easily seen from A(g, ge aX ) = 
a 2 B(X-p,X-p) mod a 3 . Note also that, since A(hgi,hg2) = A (01,(72) = A(gi/ci, (72^2) for 
/i G G and k\,k2 G K, the replacement ([16]) preserves the global and local symmetries of 
the continuum theory. 

By its origin from B [[g~ l dg)-p 1 (g~ 1 dg)-p) with d differentiating in time (space) di- 
rection, the summand with ao in the denominator (numerator) is naturally identified as 
a "kinetic energy" ("potential energy"). We now wish to extract the logarithm of the 
transfer matrix in the Hamiltonian limit clq/cli — ► 0. To do so, we must first perform a 
standard transformation, turning kinetic energy into a second-order differential operator 
while inverting the coefficient a±/ao. A crucial simplification gained from making the sys- 
tem very anisotropic is that we can ignore the potential energy terms while doing this 
transformation. (Both potential energy and the differential operator that will eventually 
represent kinetic energy carry a factor of ao, so their commutator is of negligible order a^.) 
Thus, concentrating on a single site n, setting g n = g and (7 n +e = h, and defining e by 
do/ai = a xx e, what we need to do is to work out the integral 

(rf)(n(h)) := / dg K exp (--A(h,g)) f(n(g)), (17) 
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which transfers some function /(tt to the neighboring site n + eo in time direction. 
This is prepared by the following little calculation: 



(rf)(n(h)) = [ dg K e~ A ^ ^ f (n (hg)) = [ dX e -B(*.*)A f(n(he x )). (18) 

J G/K Jv 

The first equality sign follows from the invariance of dgx under the substitution g i— > 
and from the symmetry property A(h,hg) = A(l,g) with 1 being the unit element in 
G. For the second equality sign we set g = e x , parametrizing G/K by its tangent 
space V at the origin tt(1). In the small-e limit under consideration, the integral re- 
ceives contributions only from the immediate vicinity of 7r(l). We may therefore replace 
A(l,e x ) = STre x Ae- x A/4 = STre 2X /4 = STr cosh(2X)/4 by B(X,X) = STrX 2 /2 and 
dgx by dX , the Euclidean integration measure on V . 

It remains to convert the last integral in ( |T8D into a differential operator acting on /. 
To that end, we express X in some basis {e^}j = i ) ... ) dim'P as X = ^x*ei, so that 

B(X,X) = J2x i g ij x j (19) 

with gij = B(ej, ej). By Taylor expanding / and doing the resulting Gaussian integral, we 
easily get 

(rf)(7r(h)) = (f + 4eCf)(7r(h)) mod e 2 (20) 
where £ is the second-order differential operator defined by 

(Cf)(n(g)) = \(d/dx)f(n(ge W j:x*e l ))\ x=o with X(d/dx) = J^gV^L-; . (21) 

ij 

(The indices of the metric tensor are raised by Ylij9 % ^9ji = $ l i as usual.) Note that, in 
view of its derivation from a manifestly gauge-invariant expression, C is independent of 
the choice of gauge and thus makes sense as an operator acting on functions on G/K. It 
is not difficult to see directly from (|2l]) how this gauge invariance comes about. A change 
of representative g i— > gk (k G K) amounts to the same as a coordinate transformation 
^2x l ei I— > J2y l ei = k(^2 x % ei)k~ x . Such a transformation leaves B(X, X), Eq. (|l~9|) , in- 
variant and the process of raising the indices of the metric tensor confers this invariance 
upon A. 



We can now read off the quantum Hamiltonian H from Eqs. (|13|) , flTq), ( |17|) and ( |20|) : 

a Nl ( 4 \ 
a H := -InT = — V" ( C t + a xx A(g h gi+ x ) J mod (a /ai) 2 . (22) 

The first term is site-diagonal, and for a single site is defined by the formula (|21|). The 
second term couples neighboring sites and acts on "wave functions" /(tt(^i), ir{gNi)) by 
ordinary multiplication. 

We should like to point out that the operator C is nothing but the Laplace operator 
V 2 = sdiv grad (gradient followed by superdivergence) for G/K. To see this we note that, 
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since a left translation L g : n(h) t— > n(gh) acts on the left while £ acts on the right, L g and 



C commute. Now, by a basic result of the theory of homogeneous spaces G /K [ 4(| , any 
second-order differential operator with the property of commuting with left translations is a 
multiple of the Laplace operator. Hence £ = cV 2 (c € R). The constant of proportionality 
c equals unity since A is the correctly normalized differential operator associated with 
the quadratic form defining the metric. Choosing some system of super-coordinates in 
which the metric B ((g~ 1 dg)-p, (g~ 1 dg)j>) is expressed by Yl d£ l <7ij(£)d£ J , we have for V 2 
the familiar coordinate expression 

V 2 = J(0- 1/2 E(-) |l| ^ J ^)^) 1/2 |j ( 23 ) 

ij 

where J(£) = SDet(<7ij(£)) and SDet is the super determinant. (The sign factor involving 
the superparity \i\ of £ l is due to supersymmetry.) The coordinate expression (|3|) illus- 
trates the meaning of the less familiar formula (^Tj). For practical calculations, it is much 
less useful than (pi]) and we will actually never use it. 

What can we say about the Hamiltonian H, Eq. (p22~D? The normalization of the func- 
tional integral Z = 1 translates into trace exp(—/3i7) = 1 for all f3, a xx and N\. We there- 
fore expect H to have a unique singlet ground state with zero energy and excited states 
organized into supermultiplets consisting of an equal number of bosonic and fermionic 
states, whose contributions to the partition sum exactly cancel each other. Apart from 
these special features due to supersymmetry, we expect the physics of H to be similar, 
at least qualitatively, to that of the quantum Hamiltonian of a conventional (i.e. non- 
super) nonlinear a model such as the 0(3) one. The nature of the ground state will be 
determined by the battle between kinetic energy and interaction energy (formerly called 
potential energy). The first of the two makes the field diffuse on G/K, thereby caus- 
ing strong fluctuations in the ground state when a xx is small. The interaction energy, in 
contrast, has the opposite tendency of aligning neighboring field variables. For a space di- 
mension larger than the lower critical dimension of 2, either one of the two tendencies may 
win. If a xx is small enough, the quantum fluctuations caused by Ci disorder the ground 
state. As a xx increases beyond some critical value cr* x , the interaction energy takes over, 
driving a phase transition to an ordered ground state with spontaneously broken symme- 
try. This phase transition is the Anderson (metal-insulator) transition in noninteracting 
disordered electron systems with time reversal invariance broken by magnetic impurities 
or a weak magnetic field. In d = 1 + 1, the case at hand, renormalization group analysis 
of the continuum theory ( |13"D leads us to expect that quantum fluctuations always disor- 
der the ground state on sufficiently large length scales, no matter how big is the value of 
a xx . In other words, the Hamiltonian (E21) is expected to have a nondegenerate and disor- 
dered ground state with massive excitations for any a xx . This scenario must be changed 
drastically by the inclusion of the topological term. 

3.4 Topological density and Berry phase 

In Sect. 3.3, we put the Euclidean field theory (|I3p for a xy = on a rectangular lattice 
with lattice constants ao and a± , and then took the extreme anisotropic limit ao/a± — > to 
derive the expression ( p2[) for the quantum Hamiltonian H. Our present goal is to extend 
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this procedure to the case of half-integral a xy . We can imagine several ways of discretizing 
the topological density. Recalling the Q-field expression L top = e^ u STr Qd^Qd u Q/8i we 
might be led to the scheme 

/ d 2 x STr Qd^Qd v Q — > V STr Q n (Q n+eo - Q n ) (Q n+ei - Q n ) . (24) 
Jm 

Note the conspicuous absence of the lattice constants ao and oi, indicating the topological 



(i.e. metric free) nature of the term. Although the naive choice ( |241 ) pays no attention 
to the geometric meaning of the topological density, it should work just fine when field 
configurations are slowly varying. For our purposes, this would actually be good enough 
since we will be satisfied with a quantum Hamiltonian that gives the correct low-energy 
physics near criticality. However, since triples of field variables Q separated by one unit 
in both space and time direction are coupled together, the discretized topological density 
( P1"D is hard to make any analytic progress with. Fortunately, a discretization scheme much 
better than ( p4|) exists. To formulate this scheme and understand its true meaning, we 
have to work through some differential geometry. The crucial step will be to take a curve 
S 1 — ► G/K, t i — > 7(£) and lift it to a curve t i— > g(t) on G. In this way we will be able to 



relate the topological density to the geometric phase of Berry [^2| . Discretization and the 
Hamiltonian limit will then be almost immediate. 

Before delving into our differential geometric project, let me inject a few words of 
explanation concerning the use of terminology. What we are dealing with is a supermanifold 
(more precisely, a supermanifold of the Berezin-Kostant type [[13]]). Such an object is not 
made from points - anticommuting c-numbers, unlike real numbers, can't be considered as 
points on a line - so the reader may wonder just what is the meaning of the word "curve" in 
the present context?! The answer is this. When constructing an ordinary (i.e. nonsuper) 
differential manifold one has two options: one may describe the manifold either by its 
points (i.e. as a topological space with a differentiable structure) or by its functions. Since 
a function / can be applied to a point x to produce the number f(x), functions are dual to 
the elements of the space they are functions on. Because of this relationship, an ordinary 
manifold can be reconstructed from the knowledge of its functions. In the super case the 



point option is nonexistent, which is of course why Berezin |^3] defines a supermanifold by 
its (super-)functions. This is not as horrible as it may sound, as the vast space of functions 
- on a supermanifold just as on an ordinary manifold - can be generated from a basic set 
of coordinate functions. The message following from this is that, on a pedantic level, we 
ought to be using coordinate language throughout. For example, when we talk about a 
"curve" on a supermanifold, what we should really do is to describe the object "curve" in 
terms of some set of superfunctions depending on a parameter t and taking values 
in the underlying parameter Grassmann algebra. (In the ordinary case, the £*(t) would 
result from evaluation of a set of coordinate functions x 1 : G/K — > R along the curve 
t i— > 7(t).) Clearly, our notation will be more transparent, and our general train of thought 
much easier to follow, if we say "curve" rather than "a set of functions depending on a 
parameter and taking values in ..." This is why we shall stick to point manifold language 
for the most part, keeping however in mind that everything we do should be expressible 
(and ultimately must be expressed) in terms of supercoordinates. 
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After this preamble, we consider for some smooth field configuration the integral of the 
topological density L top over some connected and simply connected region T of space-time 
M: 

</?r = J d 2 x L top . (25) 

We wish to elucidate the meaning of <pr as the geometric phase acquired by a "G/K- 
superparticle" when transported once around the curve dT — > %(g(dT)). We begin by 
writing L top , Eq. QT2]), in the form 

L top (g,dg) = ±e!»'STrA{g- 1 dpg) v {g- 1 a v g) v . 

Our first step will be to convert the integral over V into a line integral along the boundary 
dT. Since L top is induced by the two-form 

u = % - STr A{g- X dg) v A {g-'dg)^ (26) 

we will be able to apply Stokes' theorem if we can find a one-form A with exterior derivative 
cL4 = uj. We claim that if 7T+ = (1 + A)/2 and s is any section G/K — > Go, the choice 

A=-iSTnr+8- 1 ds (27) 

does the job. The argument is as follows. Because s _1 ds takes values in Lie(Go) = 
su(n, n\2n), the supertrace STr s _1 ds vanishes and we may replace the expression for A by 
A = STr As _1 ds/2z. Differentiating the relation g = s(n(g))k(g) we get s _1 ds = k(g~ 1 dg— 
k~ 1 dk)k~ 1 , which we insert into dA = % STr A(s _1 ds) A (s _1 ds)/2. Evaluation of the two- 
form d.4. on a pair of tangent vectors sX and sY gives i STr A[X, Y]/2 = i STr Y[A, X]/2. 
This vanishes for X G /C and Y G /C + V, so only the P-component of s _1 ds need be 
kept. Hence, d^4 = i STr A(s _1 ds)-p A (s _1 ds)-p/2, and since (s~ 1 ds)-p = k(g~ 1 dg)-pk~ 1 
and k~ x Ak = A, comparison with (|26|) shows that d^4 = u>, as claimed. 

Naive application of Stokes' theorem would now yield ifr = —i Jq dtSTrii + s~ 1 s, 
where s is the total derivative of s along the curve dT, which we parametrize by t G [0, 1). 
This, however, is not quite correct. There exists a difficulty of topological origin. For a 
fixed choice of gauge s, the gauge potential — i STr tt + s~ 1 s induced from A by the field 
configuration along the loop dT, may not have a smooth extension to the interior of the 
loop. This is because in the presence of topological excitations the field may wind once - 
or several times - around a Dirac-string type singularity of s as we move along dV. Correct 
application of Stokes' theorem would then require subtraction of contributions from the 
topological singularity. Fortunately, there exists a way to get around this problem, which 
is to use a gauge transformation to move the string singularities somewhere else on the 
manifold, far from where the fields are. This works locally, i.e. for a small loop V. In 
general we are forced to patch together several coordinate charts in each of which s is well- 
defined. The transition functions for going from one chart to another are given by gauge 
transformations. Thus, we are led to consider what happens if we make a t-dependent 
gauge transformation s(7(£)) i— > s(^(t))u(t). Such a transformation changes the expression 
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for (p-p by the addition of the integral — i J Q dt STr n + u(t) 1 u(t) = —i dt-fj^ STr n + lnu(t). 
This integral is quantized in integer multiples o/2-7r, by the following argument. In a ma- 
trix notation where A = Ibf ® (03) ar ® 1e is represented by ( 2 ™ ) and w(t) by 

( ^ £>(t))' We ^ ave STr 7r + lntt(t) = lnSDetA(t). The curve t ^ <y(t) traced out on 
G/K by the a model field along dT is closed, so u(l) = u(0) and A(l) = A(0). Integrating 
the total derivative we get 

r 1 d 

-i / dt— STr7r + lnw(t) =ilnSDetA(0) - i InSDet A(l), 
Jo dt 

which proves the statement by the single- valuedness of the function U(n\n) — > U(l), A 1— > 
SDet A. We conclude that the correct formula for <p-p is 



(p r = -i [ dt STr 7i + s~ 1 s + 2nq 1 (28) 
Jo 

with q some integer. The gauge- ambiguous integral § dt STr 7r_|_s _1 s is a one-dimensional 
version of what is called a Wess-Zumino terra in field theory. 

In the next step we perform one further manipulation, which will eventually allow us 
to completely integrate the expression for ip-p. (This is not yet possible in the present form 
since s and A do not commute.) Given the curve t 1— > 7(2) on the base space G/K, we 
define the lifted curve t 1— >■ g(t) on the total space G by the equations 

(A) n(g(t)) = j(t), 

(B) (g{t)- 1 m)K = 0. 



Eq. (B) has the following geometric interpretation. The time derivative g(t) is tangent to 
G at g(t). The left translated vector g{t)~ l g{t) lies in the tangent space of G at the group 
unit. It can therefore be regarded as an element of Q = Lie(G). Now recall the orthogonal 
decomposition Q = /C + V where /C = Lie(K). g(t) is called vertical if (?(t) _1 <7(£) lies in K. 
(g(t) is then tangent to a fibre of the K-bundle G — > G/K.) In contrast, Eq. (B) requires 
<7(£) _1 <7(£) to be orthogonal to /C for all t. This is expressed by saying that a solution of 
Eq. (B) is a curve on G which is everywhere horizontal (Figure 2). 

To solve Eqs. (A) and (B), we use the ansatz g{t) = s(^(t))u(t) with unknown u{t) = 
k(g(t)) G K. (A) is then satisfied automatically while (B) turns into 

(B') ^Ht)- 1 = -( S - 1 i)( 7 (t))^, 

which is a first-order differential equation determining u(t) uniquely once w(0) has been 
specified. Using (B ; ) in Eq. fl2"8p we can do the integral for ifr- 

(pr = i [ dt STr Tr+u^uit)- 1 + 2ivq = z STY tt + In %(0)) -1 %(l)) + 2ttq (29) 
Jo 
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since n + and u{t) G K commute. 

We conclude this subsection with a more "global" view of what we have done. Given 
a smooth field configuration and some value for a xy , the topological density assigns to a 
region r C M. the topological phase pr = exp(—ia xy L d 2 x L top ). The field configuration 
maps the boundary of V onto a closed curve on G/K which we parametrize by t G [0, 1) i— > 
■j(t). This curve is lifted to a curve dT — > G (not closed in general) by solving (g~ 1 g)K = 0. 
Making some choice of gauge fc(^) = s(ir(g))~ 1 g (more precisely, patching together local 
charts) we get pr — exp (a xy STr 7T+ ln/c((7(0)) _1 /c((7(l))). This way of writing pr identifies 
it as a Berry phase. The expression for pr is gauge-invariant if a xy = m G Z , in which 
case we can write 

p r = ^(M^(0))- 1 M^(l))) (30) 

where 

fi(k) = exp (m STr 7r + In fc) (31) 

is globally well-defined as a representation : K — > U(l). By the analogy with the Dirac 
monopole problem that results on setting G/K = SU(2)/U(1), we call the integer m the 
monopole charge. The calculation we have done in this subsection is an application of the 
beautiful theory |14| of connexions on the hermitean line bundles that are associated with 



a principal fibre bundle. We refrain from elaborating on this wider context here. 
3.5 Quantum Hamiltonian for o xy ^ 



The functional integral ([13]) carries the topological phase factor pm = exp(— ia xy J M d 2 x 
L top ). As was explained in the preceding subsection, we can use Stokes' theorem to convert 
J M d 2 xL top into a gauge-ambiguous line integral along the boundary of A4. The gauge 
ambiguity is unobservable if and only if a xy is an integer. We will derive the quantum 
Hamiltonian of the a model for this special case first. The generalization to the more 
interesting case of half-integral a xy will be presented afterwards. 

Let a xy = m G Z ■ For a configuration space M. without boundary, Eq. ( p8|) leads 
to Pm = exp(— irrupM) = 1? which is just another way of stating the quantization of the 
topological term. Hence, the quantum Hamiltonian H of the a model is given by Eq. (|22|), 
just as for the case of vanishing a xy . As was said earlier, H is expected to have a disordered 
ground state with an excitation gap, corresponding to localization of all states. A less dull 
situation arises when Ai does have a boundary dM. ^ 0. For the Corbino disk geometry, 
dM. consists of two pieces, an inner and an outer edge, both of which are circles ~ S 1 with 
radii R± and R<i respectively. The orientations they inherit from A4 are opposite to each 
other. Using Eq. fl28|) we get a sum of two Wess-Zumino terms, one for each of the two 
edges: 

/ d 2 x L top = —i f dr STr 7r + {s~^d T s N - s~{ x d T s^) mod27rZ (32) 
J M J 

where s\ (sjy) is s evaluated at x = R\ (x = R2). The appearance of these terms changes 
the a model dynamics at the edges only. More precisely, the time evolution of the edge de- 
grees of freedom is governed no longer by the free-particle Hamiltonian — C but by another 
operator which we are now going to construct. As before, we discretize the field theory (|T3| ) 
on a rectangular lattice with lattice vectors eo, ei and lattice constants ao, a± in time and 
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space direction, respectively. With all the preparations that were made in Sect. 3.3, we can 
focus attention on the kinetic energy in combination with the topological phase. To be con- 
crete, we consider the outer edge and take its orientation to follow the positive direction of 
time. Visiting a set of lattice sites n — »■ n + eo — > n + 2eo — > ... we encounter a sequence of 
field variables which we parametrize by 7r(g n ) — > 7r(g n e Xl ) — > n(g ri e Xl e x ' 2 ) — > ... (X{ G V). 
In the time-continuum limit ao —>■ this discrete sequence condenses into a continuous 
curve S 1 — > G/K which is closed by periodic boundary conditions for time. To get a dis- 
crete approximation to the line integral fl2~9p , we need to lift the discrete sequence of field 
variables to the total space G. Since t g(t) = g(0)e tx satisfies {g~ x g)K = for X e V, 
the lifted sequence is g n — > g n e Xl — > ff n e Xl e X2 etc. By formula (|30|) , the total topological 
phase factor accumulated along the outer edge is t L (k(g n )~ 1 k(g n Y[ i e Xi )) (time-ordered 
product). Using the representation property n(ko~ 1 k2) = ii(kQ~ k\)ii(ki~ k%) we see that 
the kinetic plus topological part of the transfer matrix r for a single site and a single step 
in negative time direction is given by 

(rf)(n(g))= [ dX e-^ x ^ l i{ki B )- x k{ 9 ^))S{^)), 
Jr 



compare Eq. (|18]). Reversal of the time direction is equivalent to the substitution \i — > pT . 
Taking the time-continuum limit e ~ ao — > we obtain for the a model degrees of freedom 
at the outer edge the single-site Hamiltonian — 4£ m /ai<7 xx where C m is the second-order 
differential operator defined by 

(C m f)(Tv(g)) = X(d/dx)^k(g)k(gexpJ2^r 1 )f^(9ew^e l )) . (33) 

x=0 



By setting [i = expmSTr7r + In = 1 we retrieve Eq. , defining the Laplace operator C 
In view of the relation between fx and the gauge connection A, Eq. (|2"T|), it should not be 
surprising that C m can be shown to be identical to the operator (V — imA) 2 (symbolic 
notation). (More precisely, if A = ^2 Aj(£)d£i, C m is the operator that results from (|23| ) 
by making the replacement d/d^ — > d / — imAj (^) . We will never use this equivalence.) 
Eq. (]3^) expresses C m through derivatives acting on the right, but for some purposes it is 
better to switch to derivatives acting on the left. This is done in Appendix A, where we 
derive the alternative formula 

(£ m f)(n(g)) = X g (d/dx)iu(k(g)k(e^g)- 1 )f(7r(e^^g)) . (34) 

x=0 

Here Xg(d/dx) = V dirn ^ g l i d 2 / dx^ dx l is defined by inverting the quadratic form B on the 
Lie algebra Q with basis {e^}. We will see in Sect. 4.6 that the single-site Hamiltonian 
— 4L m ja\o xx describes uni-directional transport through a number \m\ of edge channels. 

To get the site-diagonal part of the Hamiltonian for the degrees of freedom at the inner 
edge of the Corbino disk, we simply reverse the sign of the monopole charge m. (Recall 
that there is a relative minus sign in Eq. (|32] ) from opposite orientation of the two edges.) 
We will construct the ground states of — C m for all m in Sect. 4. 

Let us consider half-integral values of a xy next. In this case, the gauge ambiguity of the 
Wess-Zumino integral (|28|) prevents us from wiping the interior of M. clean of topological 
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phase. Decomposing the topological density a xy L top into its integer and fractional parts, 
we can convert the former part into a boundary term in the same fashion as before, but 
the latter part cannot be removed. We must treat it differently. This we do by modifying 
the continuum field theory ([13[) as follows. We organize the space-time manifold M. into a 
set of equidistant strips. The strips have a width of one lattice spacing a± and are parallel 
to the time direction. Every other strip is shaded, see Figure 3. We set the fractional part 
(= 1/2) of the topological coupling a xy to zero in the empty regions while raising it to twice 
its original value in the shaded ones. We are allowed to make such a drastic modification 
since, once again, all we are concerned with are the universal low-energy properties. These 
are determined by slowly varying fields, for which the staggered value of the modified 
topological coupling balances out. Put in formulas, we replace 0.5 J M d 2 x L top in ( |13|) by 
f M d 2 x L top where M. 1 is the union of the shaded regions. By using Stokes' theorem for 
M.\ we get a sum of line integrals: 

\ [ d 2 x L top — > -iJ^i-Y i drSTiTi+s^drSi . (35) 

1 JM I J 

The sign of the line integrals alternates because the orientation of the loops forming the 
boundary of M.\ does. 

The sceptical reader might have the impression that by modifying the topological 
coupling we are making a novel and unnecessary approximation. This impression is wrong. 



A close look at the literature (see [45[ for a very readable account) reveals that the same 
approximation is always made in this context and, as a matter of fact, there is no way 
of avoiding it. Our presentation is simply inverted compared to most published ones, 
which start from the sum of line integrals (|35|) (derived by using coherent states on an 
antiferromagnetic quantum spin chain) and end up with the topological term of the a 
model. The argument usually given is the following. One takes the alternating sum (|35p 
and breaks translational invariance by organizing the line integrals into groups of two: 



dr STr7r + (s 2Z l ) _ i a T S2z+i - s 2l 1 d T s 2 i). 



By Taylor expanding up to linear order in the lattice constant a±, one gets i J2i a i § dr STr 
7v + (d x s~ 1 d T s)2i- Because periodic boundary conditions in time direction are used, this 
expression is not changed by subtraction of a term with d x <-> <9 T , whence it becomes 



^ ai jdr STr h[ S - x d T s, s~ x d x s) 
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The expression i STr A[s~ 1 d T s, s~ 1 d x s] is recognized as twice the topological density L t0 p. 
Note in particular that the sum runs over even lattice sites only (i.e. the topological 
coupling is staggered). By replacing X] even i ~ ¥ Sail i an d taking the naive continuum 
limit, one recovers the topological term j M dxdr L top /2. 

We are finally ready to write down the quantum Hamiltonian H for a xy G Z+l/2. 
Putting everything together, we get the formula 

a x H = --^- ( C t+ C T ) +<?xxY< A (9h9l+i)- (36) 



^even I odd I J all I 
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where £ ± := C^ 1 is defined by Eq. (^3|). This expression for H replaces the formula 
and is correct modulo terms of higher order in the small parameter ao/ai- The sign of the 
monopole charge m = ±1 alternates on even and odd sites because the sign multiplying 
the line integrals (|35|) does. Eq. ( p6|) applies to case of a two-torus M. = S 1 x S 1 . For the 



disk geometry it must be modified at the edges to account for the integer part of a xy in 
the manner described above. 

So far, we have explained how to pass to the Hamiltonian formulation of the a model 
when a xy G Z +1/2. This is not yet entirely satisfactory since, in order to compute critical 
indices such as the correlation length exponent u, we may wish to study the a model 



xy 



slightly off these points. We can do so as follows |T8[. Suppose a xy = (1 — e)/2 with 
|e| <C 1. For definiteness, let e > 0. If we choose M.\, the union of the shaded regions 
in Figure 3, to satisfy the condition area(A^i) = area(.M)/2 as before, the alternating 
sum of line integrals ( p5|) will produce an excess topological phase. We can eliminate 
the excess by making area(TVli) smaller by a factor of 1 — e. This modification amounts 
to choosing an alternating sequence of lattice constants (1 — e)a\ (width of the shaded 
regions) and (1 + e)a± (width of the empty regions). It is clear that the alternating spatial 
distance between sites causes a staggering of the site-site interaction. Thus, to perturb the 
a model off the point a* y = 1/2, we may simply add to the quantum Hamiltonian (|36|) an 
alternating term r ) ^(9U 9l+i) where r ~ a xy — a 

4 Strong-Coupling Ground States 

Low-energy properties of quantum chromo dynamics in 3 + 1 dimensions, such as confine- 
ment of quarks and gluons, are beyond the reach of weak-coupling perturbation theory. 
In the late seventies, it was hoped that the incorporation of topological excitations such 
as instantons would lead to progress. These hopes, however, were frustrated by severe 
infrared divergences which make it difficult if not impossible to formulate a sensible ex- 
pansion around the dilute instanton gas limit. On the other hand, confinement occurs 
very naturally in lattice QCD, by the area law for Wilson loops in the strong-coupling 
limit. The situation is similar for the field theory ([13|) of the integer quantum Hall effect. 



Although weak-coupling perturbation theory augmented by instanton calculus [TI,r2| has 
yielded some qualitative insights, it has not led to any conclusive or quantitative statements 
about low-energy (or large-scale) phenomena such as quantization of the Hall conductance 
and the derealization transition between Hall plateaus. In this case, too, an approach 
starting from the opposite limit of strong coupling looks more promising. Diagrammatic 
perturbation theory predicts that the renormalization group flow takes a xx to small values. 
Therefore, an expansion around the sma\l-a xx (or strong-coupling) limit is more likely to 
capture the essential long wave length physics. 

The quantum Hamiltonian H, Eq. (^6|), is a sum of two terms: a\H = AHq/(j xx + 
a xx Hi. For a xx —> 0, the term H = - J2i(^tl +J ^2l+l) dominates and Hi = ^ A(g h g t+1 ) 
can be treated as a perturbation. The first step towards making a strong-coupling expan- 
sion is to diagonalize H . Since H does not couple fields on neighboring lattice sites, 
diagonalization of Hq reduces to a one-site problem. We shall denote the one-site Hamil- 
tonian by 7i s = —£ m = —(V — A) 2 , concealing the dependence on the monopole charge 
m = ±1 but making the dependence on the choice of gauge s by A = m STr 7r + s _1 ds 
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explicit. The purpose of the current section is to construct the ground state of H s for all 
m G Z • (Doing the calculation for a general value of m will enable us to discuss edge 
state transport in the Hall plateau regions in Sect. 4.6.) We will see that this ground 
state is infinitely degenerate unless m = in which case \i = exp m STr 7r_|_ In is trivial. 
States in the degenerate ground state (as in any eigenspace) of TL S transform according 
to some representation of the symmetry group G. It turns out that for m > (m < 0), 
this representation possesses a lowest (highest) weight. Our strategy will make use of this 
fact: we will first identify the lowest (highest) weight and then generate the entire space of 
ground states by symmetry transformations. Part of what we do takes place at the edge of 
current mathematical understanding. This is why we need to begin with some preparation 
in Sect. 4.1. 



4.1 Principal remarks 

We begin by recalling some of the obvious structures at hand. The basis for our consider- 
ations is the super coset space G/K = U(n, n\2n) /XJ (n\n) x U(ti|ti). G acts on G/K by 
left translation x = tv{K) \— > g • x = n(gh). This induces an action of G on functions on 
G/K by 

(L g f)(x) = f(g- 1 -x). (37) 



(Since the construction of L g f from / involves Grassmann-analytic continuation [^3[ of / 
w.r.t. the Grassmann parameters of g, f must have sufficiently good analyticity properties 



in order for L g f to be well-defined.) The action (|37| ) has the representation property 
L g h = L g Lh- 

Integration of functions on G/K is defined in the sense of Berezin |[43|| . (Again, since 
the definition of the Berezin integral involves differentiation w.r.t. the fermionic coordi- 
nates, functions must be analytic.) G/K carries with it a Berezin integration measure dx 
which is invariant under left translations, dx is unique within multiplication by a constant 
and is alternatively denoted by dgx- By the invariance of dx, 

/ f(x)dx= (L g f)(x)dx (geG). 
Jg/k Jg/~k 

Given two analytic functions fi and /2, we define their scalar product (/i, $2) by 

(/i,/ 2 )= / hix)f 2 (x)dx (38) 

J G/K 

whenever the right-hand side exists. The bar denotes complex conjugation. The scalar 
product ( [38D is preserved by left translation, 



(L g f 1 ,L g f 2 ) = (f 1 J 2 ) (geG). 

We will see later that the operator 7i s = — (V — A) 2 is hermitean with respect to it. In 
this respect, the scalar product ([JSp is similar to the scalar product on the space of square- 
integrable functions in standard quantum mechanics. There is, however, one important 
difference: the scalar product ( |3~8"1) is not positive definite and does not define a Hilbert 
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space. (See Eq. (69) below.) From a rigorous mathematical point of view, this raises 
a number of interesting questions: convergence of sequences of functions, completeness 
of the function space etc. These are the standard questions of functional analysis that 
physicists never need to worry about because they always deal with Hilbert spaces where 
these questions have been answered by mathematicians long ago. It seems to me that the 
proper definition of norm and topology for function spaces on supermanifolds is not well 
understood at present. (In any case, / don't understand it very well.) 

The loss of positivity is worrysome not just on a rigorous level but may sometimes 
cause severe problems for the practical user. As an example we mention the difficulties that 
one encounters when trying to establish a workable theory of Fourier analysis on compact 
supergroups [§6| . (By Fourier analysis we mean expansion in eigenfunctions of the Laplace 
operators.) In this case the lack of a positive definite invariant hermitean form causes 
Laplace-Casimir operators to be generically non-diagonalizablel Fortunately, the situation 
for Riemannian supersymmetric spaces G/K - as opposed to supergroups G - appears to 
be much better, in spite of the lack of positivity. The first proof of a Fourier expansion 
theorem, for the simple case of a rank-one supersymmetric space, was published in [47]]. 
The method of f|7j was subsequently extended J|8| to all of Efetov's spaces [^2|, i.e. the 
supersymmetric spaces G/K that appear in the theory of noninteracting disordered elec- 
tron systems. It was found that the problem of non-diagonalizability of Laplace operators 
is absent for these spaces. It turns out that the method invented in [^7| can be further 
generalized to prove a spectral expansion theorem for the operator 7i s . (Note that the 
hermitecity of TL S implies orthogonality of its eigenspaces w.r.t. the scalar product (|38|).) 
In the present paper we are concerned with the strong-coupling limit of the Hamiltonian 
if, Eq. ([36|). For that, all we need are the ground states of 7i s and its conjugate. We will 
therefore concentrate on these ground states and relegate the full spectral expansion of H s 
to a future complete account of my work on supersymmetric Fourier analysis. 

With superanalytic theory being incompletely developed at present, even the modest 
goal of constructing the ground state of 7i s is not entirely straightforward. Suppose we 
have found some eigenfunction of 7i s . How do we know it occurs in the spectral expansion? 
(Note that invariant differential operators on noncompact spaces are bounded neither from 
above nor from below. They have a large number of eigenfunctions that do not occur in 
their spectral expansion. Examples of such eigenfunctions are the functions that correspond 
to the scaling fields of Sect. 2.1.) On a Hilbert space the answer is easy: we simply look to 
see whether the function is (properly or improperly) normalizable. If it is not, we reject it. 
For reasons that are not entirely understood, the same principle works in the present case, 
in spite of the lacking Hilbert space property. Although the proof of the Fourier expansion 
theorem for Efetov's spaces makes no reference to the hermitean form (|38l), it turns out 
that the eigenfunctions occurring in the Fourier expansion are precisely those with "proper" 
behavior with respect to (By this I mean that eigenfunctions belonging to a discrete 

series have finite norm in a certain sense, while those belonging to a continuous series are 
normalized to a 5-function.) This then is the heuristic principle we adopt: given some 
eigenfunction of 7i s we shall declare it to be a "physical state" (occurring in the spectral 
expansion) if it is normalizable w.r.t. 
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4.2 Induced Representation 

We again begin with some motivation. Eq. ([J7[) defines an action of the symmetry group 
G on functions on G/K. This action is not yet the "good" one for our purposes, since 
it does not commute with the Hamiltonian 7i a = — (V — A) 2 . If this is not obvious, 
consider the simple example of the motion of a charged particle in the field of a magnetic 
monopole. Such a system is invariant under rotations that fix the magnetic monopole. 
However, rotational invariance is not manifest in the Hamiltonian formulation since writing 
(V — A) 2 requires making some choice of A, which inevitably breaks spherical symmetry. 
Consequently, the symmetry transformations of 7i s are not just rotations but are rotations 
followed by a suitable gauge transformation. The same happens in the present case. Thus, 
the good action of G will combine left translations with gauge transformations. 

To construct the good action we proceed as follows. For a fixed value of the monopole 
charge m, let us denote the linear space of functions which 7i s acts on by W s . The subscript 
s does not mean inequivalence of W s and W s i but is simply a reminder that the elements 
of W s are tied to some choice of gauge A = mSTr 7r + s~ 1 ds. When s is changed, the 
functions change. Thus, W s is a space of gauge- dependent functions on G/K. Since a 
change of section s i— > sk takes A into A + d(m STr 7r_|_ In k), we expect the behavior under 
gauge transformations to be governed by the one-dimensional representation \x of the gauge 
group K which is defined by \i{k) = exp (mSTr7r_|_ Ink). To confirm this expectation we 



use an idea that has been described to physicists in Refs. |41],f49|. It is to get rid of the 
gauge dependence by "working with all gauges at once" and regarding an element / G W s 
as a function F on the total space G. In detail this works as follows. Recall that a choice 
of section s gives rise to a map k : G — > K by the equation g = s(irr(g))k(g). Given k(g), 
we associate with / G W s a function F on G by 

F(g)=^k(g))- 1 f(n(g)). (39) 

Since k(gki) = k{g)k\ for k\ G K and since fi is a representation, F satisfies 

F(gh) = tiihy'Fig). (40) 

In this way, W s corresponds to the linear space, W. of functions F on G subject to the 
symmetry condition (J40Q . The correspondence is one-to-one, with the inverse of relation 
being 

f(x) = /j(k(s(x)))F(s(x)) (xGG/K). (41) 

In other words, Eqs. ( f41~D and (^9|) define a bijection ifj s : W — > W s . (An element / G W s 
is technically called a section of an associated bundle |44| .) 

Let us illustrate the above correpondence between / and F at the simple example of 
SU(2)/U(1). With the choices 

g = e i 4>°~3/2 e 'iQ°~2/'2 e i^a 3 /2 g _ e i(f}a 3 /2 e i8a 2 /2 ^(e* 00 " 3 / 2 ) = e iKa 

the eigenfunctions of the Dirac monopole Hamiltonian are D( M (« _1 ) = P^ M (0, —9, —<fi) 
= f where T> are Wigner's D-functions expressed in Euler angles 4>, 9,ip. They extend to 
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functions V J KM {g 1 ) = T>K M (—ip, —9, —(f)) = F, which are functions on SU(2) obeying the 
symmetry condition 

Vkm ((^ 3/2 )" 1 ) = e- iKa V J KM {g- x ). 

How does the monopole Hamiltonian operate on functions F G W7 Making the 
similarity transformation Ti := if)~ 1 H s if) s we get from Eq. (^3|) for 7i s = — £ m : 

(HF)(g) = -X(d/dx)F(gexp^x l e,) . (42) 

x=0 



Note that Ti preserves the symmetry condition ([I0p. In fact, assuming F to satisfy 
and using the invariance of X(d/dx) under J2 xle i l— *■ k(^2 l x % ei)k~ 1 , we get 



x=0 



(HF)(gk) = ~X(d/dx)F(gkexp^< 

= -n{k)- 1 \{d/dx)F (gexpik^e^- 1 )) = ^k)-\HF)(g). 



2 = 



Therefore, Ti is well-defined as an operator W — ► W . 

Consider now the natural representation of the symmetry group G on W by left 
translation: 

(L g F)(h)=F(g- 1 h). (43) 

(Right translations are unnatural because they interfere with the constraint (^0|).) The 
representation (f|3D commutes with the operator TC, Eq. (|4~2]). This is obvious since left 
translations act on the left while Ti, acts on the right. At this point we recognize the advan- 
tage of making the similarity transformation from 7i s to Ti: the symmetry transformations 
of Ti are simply the left translations! The good action, T g , of the symmetry group G on 
gauge-dependent functions / G W s now follows by transforming from W back to W s . We 
have T g = -0 s L £/ , i/' s _1 or, in explicit terms, 

(T g f)(n(h)) = ^{kmig^h)- 1 )}^- 1 ^). (44) 

These are the symmetry transformations of 7i s . They form a representation of G, by virtue 



of their construction from ( [43]) via similarity transformation. The action (g4|) is called an 
induced representation. It depends on \x and on the choice of gauge s. This fact is not 
made explicit in our notation. 

We claim that the action (|44|) preserves the hermitean form (|38|) : 



(T fl / 1 ,T ff / 2 ) = (/ 1 ,/ 2 ) (g G G). 

(If (•, •) were positive definite this would amount to saying that the representation is uni- 
tary.) In view of Eq. ([14[) and the invariance of dx under left translations x = tt(/i) i— > 
g~ x ■ x = 7r((7 _1 /i), our statement is obvious if we can show that 



H(k(g)) = (45) 
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holds for g E G. The relation ([±5]) is established as follows. From E g. (|3T|) and STr A = 
STr A T (superscript T denoting supertranspose) and A T = A' , we get fi(k(g)) = fx(k(g)) = 
n(k(gY). Now, by Eq. (Q) k(g) E K C G satisfies k(gy = Tjk(g)~ 1 rj and, since fi is a 
representation and fi(r]) 2 = 1 (note rj E K), fi(k(gy) = n(k(g))~ 1 , so ( f45|) is proved. 

It follows from (|45|) that the monopole Hamiltonian 7i s is hermitean w.r.t. the scalar 
product (BSD- The proof is simply a matter of writing down 



(f2,H a fi) = [ dg K f 2 (7r(g))(HJi)(7r(g)), 

J G/K 



'G/K 

using the formula (|34]) for 7i s = — £ m , making the change of integration variable g i— ► 
e -Ex ( w hich leaves (igx invariant), and then using Eq. (flop. 

4.3 Choice of gauge, and coordinate presentation 

To proceed further, we pick suitable coordinates for G/K, and we fix the gauge by choosing 
some section s : G/K — > G. Our choice of coordinates is motivated by the requirement 
that the projection tc : G — > G /K be expressed in a particularly simple form. 
An element g E G is written 

»={c d)- < 46 > 

Here, and in all that follows, the displayed matrix structure refers to advanced-retarded 
space. Each of the blocks A, B, C, D is a 2n x 2n supermatrix acting in the tensor product 
of Boson-Fermion and Extra space. They are subject to a certain set of constraints that 
follow from the defining equation for G, Eq. (^). In particular, these constraints imply 
that both A and D have an inverse in some neighborhood of the unit element of G. Since 
A = diag(lsFigi_B, — 1bf®e) in the presentation adopted, an element k E K is of the 
block-diagonal form k = * fc V Invariants for the action g i— > gk (k E K) are the 

combinations CA~ X and BD~ l . Counting shows that the number of independent variables 
in CA~ l and BD~ l equals the number of degrees of freedom on G/K. We put Z := CA~ X , 
Z := BD~ X and take Z, Z as local coordinates for the coset space G/K. Thus, f(n(g)) 
will from now on be written f(Z, Z). (For clean mathematics we should work with several 
coordinate charts. We will however get away with using just one.) Recalling the defining 
equation for G one can show with a little calculation that Z and Z are related by Z = Z't^ 
where r 3 = (a 3 ) B F ® Ie- 

The meaning of the coordinates Z, Z is best understood by looking at some examples. 
Let us specialize to n = 1, in which case G/K has a two-sphere S 2 ~ SU(2)/U(1) ~ 
U(2)/U(l) x U(l) for its FF-space. By projecting the pair of supermatrices (Z, Z) on the 
FF-sector, we get a pair of complex numbers (—z, z). (The minus sign is dictated by the 
unitarity condition Z = ZV3.) The special element A is replaced by the Pauli matrix (73. 
Parametrizing SU(2) by Euler angles 

g = exp(z</>a 3 /2) exp(z#a 2 /2) exp(#<7 3 /2) = ("p & ) > ( 47 ) 
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we get z = (3/a = tan(0/2) exp(— i<j)). It follows that, if n = (m, 712, 713) = (sin 9 cos 4>, sin 9 
sin 0, cos 0) is a unit vector parametrizing the two-sphere S 2 , z has the expression z = 
{n\ — in<i) /(l + n 3 ). In other words, z is the usual complex coordinate associated with the 
stereographic mapping of S 2 ~ SU(2)/U(1) onto the complex plane C • Consider next the 
noncompact analog space U(l, 1)/U(1) x U(l) ~ SU(1, 1)/U(1), which is the BB-space of 
G/K for n = 1 and is isomorphic to the two-hyperboloid (or hyperbolic plane) H 2 . For this 
we replace (Z, Z) — ► (z, z) (without the minus sign). Changing to Euler angles for SU(1, 1) 
by making in the exponent of the middle factor in Eq. ( fj7|) the substitution i&2 — > ci, 
we get z = tanh((9/2) exp(— i(p). H 2 can be parametrized by a hyperbolic unit vector 
n = (ni, 712,713) = (sinh6> cos cf>, sinh6> sin</>, cosh 9). As before, z = (711 — iri2)/(l + 713). 
What this means is that z is the complex coordinate of the Poincare disk model (\z\ 2 < 1) 
of H 2 1 150|| . Finally, we mention that Z and Z coincide with the matrices W21 and — W12 
used for the parametrization of Q in Ref. [12] for n = 1. To verify this, one expresses 
Q = gAg~ x in terms of Z, Z and compares the result with Eq. (3.2) of that reference. 

Let us now fix s : G/K — > Go = SU(n, n\2n). For various reasons that will be- 
come clear as we go along, the optimal choice is s(Tr(g)) = (gAg~ 1 A) 1 / 2 with coordinate 
expression 

!)( (1 -f rl/2 m> 

It is easy to verify that (i) s(Z, ZV3) G Go and (ii) 7r o s = id. (These properties are 
sufficient in order for s to be a section of the bundle G — > G/K.) Specializing to the 
FF-sector for n = 1 and using the Euler angle parametrization (|1"7|), we see that the choice 
s(7r(<7)) = ((7C7" 3 (7 _1 a3) 1 / 2 corresponds to the second of the two sections in ([15]). 

We now wish to re-express the central formulas (|33|) and (|44]) in terms of the complex 
coordinates Z and Z, starting with the latter. For that we need to figure out how left 

translations 7t(h) 1— > ix(gh) act on Z and Z. We put h = f , ^, J and characterize the 

coset 7r(/i) by Z = C'A' 1 and Z = B'D' 1 . For the product g/i with g given by Eq. (|4"6| ) 

we write gh = ( „ D „ J. The complex coordinates of the left translated coset ir(gh) are 

denoted by g ■ Z = C" A" 1 and g ■ Z = B"D" 1 . Carrying out the matrix multiplication 
we find 

g-Z= (DZ + C)(A + BZ)-\ 

(49) 

g-Z= (AZ + B)(D + CZ)~ 1 . 

We mention in passing that these relations generalize the rational action of SL(2, C) on C 
by z I— > ffqrf • By solving Eqs. (^9|) for Z, Z we get the inverse transformations 

■ Z = (D — ZBY X (ZA - C), 
# ■ Z = (A - ZC)-\ZD - B). 



The other ingredient on the right-hand side of Eq. fll3| ) is the multiplier built from fj, and 
k(-). Unfortunately, the derivation of its explicit expression in terms of Z, Z requires a 
certain amount of algebra. To avoid cluttering the main text with a lot of formulas, I have 
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relegated these details to Appendix B. Using the information given there, we arrive at the 
following coordinate expression for the representation ([14]): 



1/2 

(T g f){Z, Z) = (SDet g)^ 2 SDet ( ) fig- 1 ■ Z, g- 1 ■ Z). (51) 

(Eq. (fJB|) is understood.) Next, we write down the coordinate expression of the second- 
order differential operator 7i s = —C 171 defined in (|3^). Again, the calculations are some- 
what involved which is why they are done in Appendix B. The result however is quite 
transparent and can be stated as follows. Let X, X be a pair of supermatrices with 
the same dimensions and symmetries as Z, Z. If we define a differential operator d x x 
analogous to X(d/dx) in fl55| ) by 

d xj} = Y(-) W ^L- , (52) 

x ' x ^ dX ll 3X ll y ' 

the operator H s has the coordinate expression 

/ ~ \ rn/2 

(H s f)(Z,Z) = -d xx SDet[ j f{(Z + X)(l + ZX)-\...) ^ y . (53) 



1 + ZX 



x=x=o 



The second argument of / on the right-hand side is obtained from the first one by Z <-> Z 
and X <r-> X. 

Since the formula fl53|) for the monopole Hamiltonian 7i s may look unfamiliar we 
turn, once again, to the prototypical example SU(2)/U(1) for illustration. In this case it is 
well-known |l(| that, for m = 1, 7i s has a spin 1/2 ground state spanned by the functions 
T>_y 2 _ 1 / 2 ( s ~ 1 ) = cos(0/2) and T>_^, 2 1 / 2 ( s ~ 1 ) = sin(0/2) exp(— itfi). To verify this from 
the present formalism, we use 

(H s f)(z,z) = -JL (iz^V^f ff^ , c.c.) , (54) 

as results from (53) by specialization to the FF-sector. (The multiplier got inverted because 
SDet puts the determinant of the FF-block in the denominator.) The expressions for the 
functions cos(#/2) and sin(0/2) exp(— i<jj) in terms of the complex coordinates z, z are 
(1 + zz)~ x l 2 and z(l + zz)~ x l 2 . From Eq. fl54] ) these are easily seen to be eigenfunctions, 
for m = 1, of 7i s with eigenvalue 1/2. More generally, for any positive monopole charge m 
the functions z k (l + zz)~ m l 2 (0 < k < m) are seen to be eigenfunctions with eigenvalue 
m/2. They form an irreducible spin m/2 multiplet of SU(2). Indeed, by specializing 
Eq. (^T[) and making the similarity transformation f(z,z) = (1 +zz)- m / 2 <p(z) we get the 
transformation law 

{ r ww = (a + ^(^). (55) 

\-p &) 
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From this we read off the representation of the generators cx_|_ = (o"i + zo"2)/2, cr_ = 
{<J\— icr 2 )/2, and (73 by the differential operators — z 2 d z +mz, d z , and 2zd z —m, respectively. 
It can be shown that the (m + l)-dimensional space of functions f(z, z) = (l + zz)~ rn / 2 (p(z) 
with holomorphic f{z) = ^^=0°^ ls the ground state of H s . (Note that in the zero 
curvature limit m — > 00 we recover the lowest Landau level wave functions f(z,z) = 
exp(— zz)(p(z) of the Hamiltonian —(V — A) 2 on the complex plane in the symmetric 
gauge A = (zdz — zdz)/2.) The representation of SU(2) on holomorphic functions has 



been used |5lj in numerical studies of the fractional quantum Hall effect, (z = v/u in the 



notation of Ref. J51[.) Eq. fl5"5]) applies to the case of positive m. When m is negative, we 
get a ground state spanned by antiholomorphic functions. 

All of the above can be transcribed without difficulty to the noncompact analog space 
SU(1, 1)/U(1). To cut things short, we just mention that the ground states z k (l + zz)~ m / 2 
(0 < k < m) translate into the ground states z k (l — zz) m l 2 (0 < k < 00). The latter 
transform according to a unitary infinite-dimensional lowest-weight representation of the 
noncompact group SU(1, 1). 

4.4 Holomorphic sections 

At the end of the preceding subsection we recalled the fact that the degenerate ground 
state of the monopole Hamiltonian 7i s for S 2 and monopole charge m > is spanned by the 
functions (1 + zz)~ rn ^ 2 (p(z) with holomorphic (f(z) = Yl J k=o c k zk - This observation holds 
the clue how to make progress with 7i s in the general case, Eq. (|53|), and nonnegative m. 
The functions (1 + zz)- m / 2 for S 2 and (1 - zz)™/ 2 for H 2 , generalize to SDet(l - ZZ) m / 2 
for our quantum Hall model space G/K. Let us therefore set v m (Z, Z) = SDet(l — ZZ) 171 / 2 
and make a similarity transformation from 7i s to 71+ = v ^'H s v rn . The transformation is 
carried out with the help of the identity 

1- (Z + X)(l + zx)-\z + X)(l + zx)- 1 = (i-zz)(i+xz)- 1 (i-xl > )(i + zl > )- 1 . 

Substituting / = VmW into Eq. ( |53|) we obtain 



in 



(1 _ xx) 1 / 2 

(H+w)(Z,Z) = -d x ^SBet[ { I ) w((Z + X)(l + ZX)-\...) (56) 

1 -|- Z X I X—X—O 

What happens for holomorphic w(Z, Z) = f(Z)7 The crucial feature of Eq. (|56|) is the 
absence of any dependence on X in the denominator of the multiplier. With this obser- 
vation, and recalling the definition of d x in (|52|), we immediately find f(Z) to be an 
eigenfunction of 7i+ with eigenvalue 

-d x x-SDetfl -XX) m / 2 . (57) 

x > x v ' x=x=o 

By the property of perfect grading, i.e. the equal number of bosonic and fermionic degrees 
of freedom, the eigenvalue (|57|) vanishes. We thus arrive at the important conclusion that 
functions / of the form 

f(Z, Z) = SDet(l - ZZ) m/2 ip(Z) 
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with holomorphic tp are zero modes of the Hamiltonian 7i s for monopole charge m > 0. For 
the general reason of stability of the functional integral (|13|) , we expect only eigenf unctions 
with nonnegative eigenvalues to contribute to the spectral expansion of 7i s . (The operator 
7i s does have a vast number of eigenf unctions with negative eigenvalues. These however 
are not normalizable and do not occur in its spectral expansion.) We thus perceive the 
central role played by holomorphic functions: they span the degenerate ground state of H s 
not only for SU(2)/U(1) but for our quantum Hall model space G/K (and other spaces 
G/K) as well. The story does not end here, since the requirement of decent behavior 
under transition to another coordinate chart imposes a constraint on <p(Z). (Recall that 
for SU(2)/U(1) the degree of the polynomial <p(z) must not exceed m.) We can avoid 
getting into the boring technicalities of switching charts by first identifying some special 
zero mode and then generating the entire space of zero modes by the action fl5"l~|) of the 
symmetry group G. This is what we will do in Sect. 4.5. To prepare the second step, we 
now specialize the representation (|5l| ) to the case of holomorphic sections. 

As is shown in Appendix B, by making the similarity transformation T+ = v^T g v m 
and keeping only the dependence on the holomorphic coordinate Z, we obtain from ( |5lD : 



(T+<p)(Z) = (SDet g) m SDet(£> - ZB)- m <p{g- x ■ Z). (58) 

(Again, Eq. ( [46]) is understood.) Eq. (|5S| ) gives rise to a representation of Q = Lie(G) by 
the usual procedure of linearization around the group unit, and this representation extends 
to the complexified Lie algebra Qc = gl(2^ 5 2n) by linearity. The motivation for going to 
the complexification is that we desire the convenience of working with the usual raising and 
lowering operators. (These do not exist inside Q.) Making the extension and then returning 
to group level, we get a representation of the complexified group Gc = GL(2n, 2n). It is 
given by the formula (|58|) without any change. (In fact, the condition g~ x = wg^r] defining 
the pseudo-unitary group G is never used in the derivation of (|58|). All we need is the 
existence of the inverse # _1 .) Having made the extension to Gc we can now write down 

the action of raising operators ( q ^ j and lowering operators ( c ° \ at the group level. 

These are seen from (^Hj) and (|5l|) to be represented by 

(T+ v?)(Z) = SDet(l-ZB)-"V((l-ZB)- 1 Z), (59a) 
U l / 

(T+ <p)(Z)=<p(Z-C). (59b) 
\c l) 

Another special case of importance is the action of elements of K<c: 

(T+ <p)(Z) = SDetA m cp(D- 1 ZA). (60) 

Vo d) 

All of the above applies to the case of positive monopole charge m > but is easily 
transcribed to m < 0. In the latter case the zero modes of TL S are functions / = V- m (p 
with antiholomorphic f(Z). Carrying out the similarity transformation T~ = vZ m T g V- m , 
we obtain 

(T-<p)(Z) = SDet(A - ZCr^g- 1 ■ Z). (61) 
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4.5 Zero modes and coherent states 



Our results of Sect. 4.4 reveal the ground state of H s for m 7^ to be vastly degenerate: 
any function / of the form / = v m (p (/ = V- m ip) with holomorphic (antiholomorphic) (p 
is a zero mode of 7i s for m > (m < 0). However, not all of these are physical states and 
contribute to the spectral expansion of TC S . It is the purpose of the current subsection to 
describe which of them do. 

We put m > 1 for definiteness and consider the special vector v m (Z, Z) = SDet(l — 
ZZ) m / 2 obtained by setting <p = 1. Our first task will be to show that v m has finite length. 
To do this with a minimum amount of effort, we put v m in coordinate- free form. Using 
g = s(7r(g))k(g) and the expression ([48]) for s, we get 

v m (7t(g)) = SDet (2 + gAg^A + AgAg- 1 ) " m/2 . (62) 



Note that the coordinate-free formula (|6^) guarantees good behavior of v m under transition 
to another coordinate chart. Aside from decaying sufficiently fast at infinity - so that 
(v m , v m ) in fact exists - v m has the special property v m (k-x) = v m (x) for k G K. Functions 
with this property are called IZ-radial. For / any such function (with the additional 
property of being integrable) , the so-called Parisi-Sourlas-Efetov-Wegner supersymmetric 
integral formula states that 

/ f(x)dx = f(o) (63) 

J G/K 

where o = 7r(l) is the unique element of G/K which is left invariant by K. (Incidentally, 
Eq. (|6^) fixes a natural normalization of dx.) A pedestrian proof of (|63|) for n = 1 can be 
found in p^] . The general proof has been worked out in [[53 1. Application of theorem (|63p 



to / = v m v m yields (v m ,v m ) = 1 since v(o) = 1 from Eq. (|62"D. This says that v m is a 
physical state. 

Given v m , we can apply a symmetry transformation T g to form the vector v g = T g v m . 
Since T g preserves the scalar product (|38| ) for g G G we have (v g ,v g ) = (v m ,v m ) = 1 
and, since T g commutes with 7i s , v g is an eigenstate of H s with zero energy. Under the 
assumption of irreducibility of the representations entering the spectral expansion of 7i s , 
we can generate the entire zero-energy sector in this way. 

For the explicit construction of all zero-energy states it is convenient to go to the com- 
plexification Gc and work with raising and lowering operators, which are norm-changing. 
A unique feature of v m is that it satisfies the equations 

1 o\Vm=v m and T k v m = fJ.{k)v m , (64) 
c 1) 

which are easily verified using (59b), (|60| ) and ([T[])- Thus, v m is a lowest-weight vector 
carrying the one- dimensional representation \i of K. To see how the raising operators act 
on v m we set (p = 1 in Eq. (59a) and find 

(2> 1 B ^Vm)(Z,Z)=SBet(l-BZ)- m v m (Z,Z). (65) 
1 
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The space of states that can be reached by repeated application of a raising operator to 
v m will be denoted by V m . It is called a lowest- weight module of G. 

A vector of the form ( |65"D , obtained by applying a raising operator exp y q q j 

parametrized by X to a lowest-weight vector, is commonly called a coherent state. To 
enhance the readability of our notation, we temporarily promote the subscript of T to 
argument and write 

Vx: = T {l X i) Vm - 

We will now calculate the coherent state overlap function J\f(X, Y) = (vx, Vy). Knowing 
it we will be able to deduce the norms and overlaps of vectors by Taylor expansion. Since 

( xi q j is an element of Lie(G) and T g for g E G preserves the scalar product, 

M(xX) = (v, n ,T(_ T [ x , °)r(; \)v m) . 

Using the representation property and carrying out the matrix multiplication, we get 

fV(X,Y) = {v m ,T^_^ % 



Y 

-r 3 xw) Vmj 



(66) 



To calculate this, we use the following trick. There exist matrices B, C, k+ and fc_ such 
that 

1 Y \ _ f 1 B\ f k+ \ / 1 0' 

-r 3 Xi 1 -r 3 X^YJ ~ \0 l)\0 k-J \C 1 

From (|6^) we have 

(v m ,T(l B )t(^ °)t(1 °)v m ) = ^(f fc n + °))(v m ,v m ). 



1 J \ k-J \C 1 J ' ' s \ k 

Hence, Af = exp (m STr In k+) = (SDet /c+) m , and solving Eq. ( |66D for k+ we obtain 

/V(X, Y) = SDet(l - T 3 X^Y)- m . (67) 

From this result in combination with Eqs. (59) and Q60D, all desired information about the 
state vector content and norms and overlaps of vectors in V m can be deduced by Taylor 
expansion. 

Let us now take a close-up look at V m for the minimal model n = 1. The base 
manifold of the super coset space G/K for n = 1 is the direct product of a two-sphere 
S 2 ~ SU(2)/U(1) and a two-hyperboloid H 2 ~ SU(1, 1)/U(1). The Cartan algebra of 
Lie(G) is four-dimensional in this case and can be taken to be generated by 

Hi = (E u )bf ® (03) ar (i = 1, 2), H 3 = (a 3 ) BF <8> 1 AR , H 4 = 1 BF ® 1 AR . 

The eigenvalues (also called weights) of the generators Hi and H2 will be used to label 
vectors in V m . Hi (H2) is the generator of the U(l) subgroup of rotations fixing the north 
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pole on H 2 (S 2 ). Introducing B^ = (Eij) BF ® (E 12 )ar (hj = 1,2), we define differential 



operators B 



V m 



V m by 



d 



&ijf = dT(Bij)f = —T exptB ..f 



t=o 



(68) 



where t is a commuting or anticommuting element of the parameter Grassmann algebra 
according to whether B t j is even or odd, respectively. Their coordinate expressions are 
easily found from Eq. (59a). Acting repeatedly with elements of the set {Bij} on v m , 
we get v m times a polynomial in the complex coordinates Z. By applying Bij once, we 
raise the degree of the polynomial by one unit. (We say that v m is the vacuum and one 
action of an operator Bij creates one quantum of excitation.) It is quite straightforward to 
construct all polynomials that are generated in this way. The degree in the anticommuting 
variables Zij (i ^ j) cannot exceed one, and the degree in the variable Z22 (coordinatizing 
the compact FF-space) turns out to be < to. Omitting the details of this straightforward 
calculation, we now give a summary of the results found. For every pair of integers p, q in 
the range p>m, —to < q <m (a pq ), p > to + 1, — m + 1 < q < m — 1 (f3 pq and 7 P9 ), and 
p > to + 2, —to + 2 < q < m — 2 (d pq ), we define the state vectors 



*P9 



(p-m)/2g(q+m)/2 



22 



V r 



Ppq — °11 °22 Bl2V m , 

Ipq - #11 #22 B 2 lV m , 

J — »(p-m-2)/2»(qr+m-2)/2/„ „ m -l R k> v, 
a pq — °11 °22 (.O12O21 — m O11O22 J^r 



These vectors form an orthogonal basis of V m . (The allowed values of p and q differ by two 
units in each case since exponents must be integer- valued.) All of them are eigenvectors 
of dT(Hi) and dT(H 2 ) with eigenvalue p and q, respectively. Furthermore, dT(H 3 )a pq = 
dT(Hs)d pq = and dT(Hs)(3 pq = +2(3 pq , dT(Hs)j pq = —2-f pq . All vectors are annihilated 
by dT^H^). The vectors a pq , (3 pq , 7 p(J , d pq for fixed p and variable q form an SU(2)-multiplet 
with spin S = to/2, (to — l)/2, (to — l)/2, (to — 2)/2, respectively. When g is held fixed 
and p is varied, they form a unitary lowest-weight representation of SU(1, 1), again with 
spin S = to/2, (to — l)/2, (to — l)/2, (to — 2)/2, respectively. The norms \f\ 2 = (/, /) are 



A pgl 



= iV 



p — to q + m 



IA 



pgi 



-|7p<j| 2 = m N 



p — to — 1 g + to — 1 



|d P9 | 2 = (to 2 - 1) iV 



2 ' 2 
p — to — 2 q + m — 2 



(69) 



where is defined by 



A(p, g) = p! 



|2 / m + p - 1 \ |2 / to 



39 



The d pq are seen to be null vectors for m = 1. The weight diagram for m = 2 is shown in 
Figure 4. 

We thus have a complete and detailed description of the lowest- weight module V m 
for m > and n = 1. The transcription to m < can be made very simply by complex 
conjugation, taking the lowest-weight vector v m into the highest-weight vector v- m and 
holomorphic functions into antiholomorphic ones. There remain two important questions 
that are left open by our analysis, i) Does the module V m really exhaust the space of 
zero-energy states? (It is conceivable in principle that Vm might be an invariant subspace 
of a bigger indecomposable representation. If this were the case, there would exist other 
zero-energy states that cannot be reached by applying raising operators to the lowest- 
weight vector v m .) ii) Are the excited states of Ti s separated from the ground state by 
a gap? Unfortunately, giving a convincing answer to these questions is beyond the scope 
of this paper. Without proof I claim that the answer to both questions is yes. As was 
explained in Sect. 4.1, I hope to publish the proof in a future complete account of my work 
on supersymmetric Fourier analysis. 

4.6 Edge dynamics and Hall conductance 

Our derivation of the quantum Hamiltonian of the nonlinear a model goes through even 
when the configuration space M. of the 2d electron gas has a boundary. The only condition 
on A4 is that one of its two directions be periodic. We will now use this flexibility for a 
semi-realistic calculation of the quantized Hall conductance in the plateau regions. The 
experimental setup we imagine is shown in Figure 5. We take for A4 a disk [i?i,i?2] x S 1 
and attach various voltage and current probes to its outer edge. 

Let us begin by summarizing what phenomenological theory M tells us about this 
kind of system. When the magnetic field is located in a Hall plateau region, all stationary 
electron states in the interior of the disk are localized. Transport proceeds via edge states 
and is effectively one-dimensional. The number M of edge channels equals the number 
of Landau levels below the Fermi energy. Edge channels are directed, i.e. electrons can 
propagate in one direction only. The direction of propagation is determined by the sign 
of the magnetic field (E x B drift). For ideal probes, the conductance coefficients g pq are 
given by 

921 = 932 = ■■■ = 9in = M, 

9i2 = 923 = ■•• = 9ni = 0, (70) 

911 = 922 = ■■■ = 9NN = ~M, 

if the ordering of the contacts follows the orientation induced by the magnetic field. All 
other conductance coefficients vanish. We will now show how the relations (|T^) follow from 
the effective field theory (p~3|) . 

The disorder averages of the g pq are given by certain correlation functions of the field 
theory whose form is specified below. To take the Hamiltonian limit of the field theory, 
we assign to the radial and angular directions of M. the roles of space and time (or, 
rather, imaginary time), respectively. By proceeding as in Sects. 3.3-5, we get a quantum 
Hamiltonian H governing the evolution in time (i.e. around the disk) of a discrete set of 
degrees of freedom. Each of these takes values in G/K and is located on one site of a 



40 



lattice in radial direction. The degrees of freedom will be referred to as G /K-superparticles 
(SPs). We are going to consider values of the magnetic field at the center of a Hall plateau 
region, where the Hall conductivity a xy is an integer m G Z • The dynamics of the SPs 
in the interior of the disk is then governed by the Hamiltonian specified in Eq. (^2|). In 
contrast, the two SPs at the inner and the outer edge sense the additional field of a fictitious 
"magnetic monopole" - the remnant of the topological density Ltop - with charge ±m. 
The sign of the monopole charge is determined by the sign of a xy and the orientation of the 
integration contour for the corresponding Wess-Zumino term § dt STr 7r+s _1 s, see Figure 
5 and Sect. 3.4. For definiteness, we take m < and m > for the inner and outer edge, 
respectively. 

Because of the influence of the probes, the SP at the outer edge gets perturbed in 
some way every time it passes a contact. (Qualitatively, passage past the contact collapses 
the wave function and puts the SP close to the origin of G/K; see Eq. (71) below.) When 
o xx is large (high Landau level), the edge SP passes the perturbation on to its neighbors by 
the action of the coupling term of the Hamiltonian fl22|) , thereby producing a complicated 
many-superparticle problem which we are unable to solve. Nevertheless, we can deal with 
such a situation by the following argument. We expect that there exists a mass gap, that is 
to say, a finite distance over which correlations decay. (This expectation comes from weak- 
coupling perturbation theory for the field theory (|13|) , as was mentioned earlier. Ultimately, 
we would like to establish the existence of a mass gap by direct calculation for the quantum 
Hamiltonian H. This however remains as a project for the future. Here we will simply 
assume the existence of a gap and see what we can deduce from it.) When the length over 
which correlations decay is much larger than the lattice spacing, we are facing the many-SP 
problem we are unable to solve. However, we can use the renormalization group to improve 
the situation. Let us imagine changing the short-distance cutoff of the field theory ( |T3|) . 
(We may think of this cutoff scale as being the lattice constant a\ for H.) As we increase 
the cutoff, the correlation length (being a length) becomes shorter when measured in units 
of the cutoff. Continuing the process, we eventually get a correlation length of order unity 
or less. In this limit, the SPs propagate independently of each other, and the perturbation 
caused by the probes is not transferred to the interior of the disk but affects only the 
SP at the outer edge. With the assistance of the renormalization group, we thus get a 
one-SP problem. This problem we can solve. (Of course, our use of the renormalization 
group is naive and overly simplistic. No argument from the renormalization group should 
be trusted without further work when the correlation length is of the order of the lattice 
spacing or less. Nevertheless, the independent edge superparticle approximation is not 
without justification, for this approximation corresponds to transport occurring only via 
edge states localized near the sample boundary, see below.) 

In the independent edge-SP approximation, the average conductance coefficients re- 
duce to the dynamic one-SP correlation functions Tr \\ i Oi exp (—(U — ti-i)H s ) with each 
operator insertion Oi representing one contact. 7i s = —C m is the one-SP Hamiltonian, 
Eq. (j33|). The imaginary time intervals are given by ti — ti-i = 4L^_i/ ' a\o xx where a xx 
is the renormalized longitudinal conductivity, and L^_i/ai is the distance between the 
(i — l) th and the i th contact, measured in units of the lattice constant a±. 

We now need the expressions for the a model operators Oi when the conductance co- 
efficient is g pq . It will be sufficient to consider the case p ^ q since the diagonal coefficients 
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g qq are determined through the off-diagonal ones by current conservation (52 p g pq = 0). 
To avoid getting into lengthy calculations that would detract from our main purpose, we 
will give a heuristic derivation. Let us first consider a diffusive quantum dot in a weak 
magnetic field (unitary universality class). In this case the formulas for the operators Oi 
can be extracted from Appendix B of Ref. []54| . Expressing them in terms of the complex 
coordinates Z, Z we find 

O p = M p Zn SDet(l - ZZ) M * (current), 
O q = M q Zn SDet(l - ZZ) M " (voltage), 
O r = SDet ( 1 - ZZ) Mr (none) , 

where Mj is the number of scattering channels in lead i. (We are assuming the ideal 
limit of unit transmission into every open channel.) From these equations we can guess 
the analogous expressions for quantum Hall systems. The main novelty that occurs for 
a strong magnetic field as compared to a weak field is spatial separation of incoming and 
outgoing channels in the leads. Figure 6 depicts the situation when the incoming and 
outgoing channels are located on the left and right edges of the lead, respectively. For an 
ideal contact there are no direct transitions between the incoming and outgoing channels 
of the same lead. As a result, every contact operator must separate into two factors, one 
each for the left (L) and the right (R) side of the contact. For example, the expression for 
a voltage contact factors as 

O q = M q (Z L ) n SDet(l - Z L Z L ) M ^ 2 x SDet(l - Z R Z R ) M ^ 2 . 

The number of edge channels in the leads equals the number of occupied Landau levels 
w = o~ xy . Making the identification M p = M q = ... = m and changing to coordinate- free 
notation we obtain 

O p = v m (BnV m , ■) (current), 

O q = Biiv m (v m , •) (voltage), (71) 
O r = v m {v m ,-) (none), 

where v m (Z, Z) = SDet(l - ZZ) m / 2 , (B n v m )(Z,Z) = mZ n SDet(l - ZZ) m / 2 , and (•, •) 
is the inner product (p8|). 

After these preparations, we are ready to proceed to the essential point of this sub- 
section. We are going to calculate the correlation function Tr Y\ i O^exp (— (£j — ti-i)TC s ), 
using the expressions (|7l|) for the operators Oi. The following two observations enable us 
to do the calculation with ease. First, since m > 0, the ground state of 7i s is spanned by 
functions / = v m (p with holomorphic (as opposed to antiholomorphic) if. The space 
of such functions is denoted by V m as before, and we set v = v m for brevity. The 
state Bnv lies in V m . Second, when inserting a complete set of intermediate states in 
the expression for the correlation function, we may restrict the sum over states to V m , 
since any contribution from the excited states of 7i s will be exponentially small of order 
exp(— const x AL/a xx ), where aL is the distance between neighboring contacts. This al- 
lows us to replace exp( — (t{ — U-i)T-C s ) by the projection operator onto V m , denoted by 
P. 
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We must now distinguish cases. Let us first consider g pq for p = q + 1. As the time 
evolution outside the interval from t q to t q +\ does no more than project onto v, we obtain 

(g q+hq ) = {Buv.e-^-^'Bnv). 

Making an exponentially small error by replacing exp (— (t q+ i — t q )H s ) by the projec- 
tor P, and using PB\\V = Buv, we get (g q +i !q ) = (Bnv,Buv). Since (B±iV,Bnv) = 
|a m+ 2,- m | 2 = m by the first of Eqs. (|69| ), we arrive at (g q+ i jq ) = m, which is precisely the 
result expected from the phenomenological picture of transmission through edge states. 
What happens when we interchange p <-> q? As before, the time evolution outside the 
interval [t q -i, tg+2] projects on v, but now (g q>q+ i) = (v, PBuv)(v, Pv)(BuV, Pv) vanishes 
by the orthogonality relation (v, B\\v) = 0. Similarly, all (g pq ) with p and q separated by 
more than one unit vanish by orthogonality. 

In the same way, we can calculate the conductance fluctuations. To get the second 
moment, what we do ||8] is to replace Z\\ (Zn) in the expression for a voltage (current) 
contact by ZnZ 2 2 {Z\\Zti)- This amounts to substituting BwB^iV for B±±v in (|7T| ). The 
factorization property (B11B22V, BwB^iV) = (i3nt>, £>iit>)(£>22t>, $22^) leads to va,r(g pq ) = 
0, so conductance coefficients do not fluctuate. All of the above results hold modulo 
exponentially small corrections and for ideal contacts and leads. 

As it stands, our argument applies only to integral values of a xy . It extends to other 
values by Khmel'nitskii's flow diagram |55[ or, more precisely, by assuming the nonlinear 
a model at (a xx , a xy ) = (0, m G Z) to be an attractive fixed point of the renormalization 
group - with the basin of attraction being the strip m — 1/2 < a xy < m + 1/2. (In the 
introduction to this paper we reviewed the argument why the critical nonlinear a model 
at a xy = 1/2 (mod 1) cannot be renormalizable but must flow under renormalization to 
another, conformal invariant field theory. There is no doubt, however, that the noncritical 
system will be attracted to the a model with a xx — > and a xy an integer.) 

In summary, for a xy close to m G N we have argued that the conductance coefficients 

92i = 932 = ... = m 

are integers, while all others (except for the diagonal ones, which are not easily accessible 
by direct calculation but are determined by current conservation) vanish. By the reasoning 
of Ref. ||, this explains the experimental observation of the integer quantum Hall effect 
for systems with ideal probes. 

5 Quantum Spin Chain 

Apart from demonstrating how quantization of the Hall conductance arises in the a model 
framework, the calculation of Sect. 4.6 serves another useful purpose by throwing much 
light on the physical meaning of the holomorphic zero-energy module V m (m > 0) and 
its antiholomorphic dual V_ m . The lesson we have learnt is this. The voltage contact q 
injects probability flux into the m edge states emanating from the lead q. In our field- 
theoretic formalism the process of injection corresponds to preparing the edge superparticle 
(SP) in the quantum state B\\v m . This state, being a zero-energy state of the one-SP 
Hamiltonian — C m , evolves in imaginary time without attenuation. A rather different 
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behavior obtains when we replace C m by its dual C~ m , which corresponds to reversing the 
sign of the magnetic field. In this case B\\v m lies outside the zero-energy space of — £ _7n , 
and therefore exp(tC~ m )Bnv m decays exponentially as t increases, with the decay time 
being determined by the smallest excitation energy of —C~ m . (Of course, here as always 
"time" is one of the two directions of the 2d quantum Hall system.) We thus see that 
evolution of edge SP states over long times (i.e. distances) is directed by the magnetic 
field. This directedness is the field-theoretic equivalent of E x B drift along a boundary. 

Having understood this, we now turn to the critical value a xy = 1/2 (mod 1) corre- 
sponding to a point of transition between Hall plateaus. The quantum Hamiltonian, H, 
of the nonlinear a model for such a point has been given in Eq. fl36|). H differs from the 
Hamiltonian for a plateau region by the alternating sequence of monopole Laplacians 
replacing the plain Laplacian C This difference is crucial: while — £ has a nondegenerate 
ground state with a gap, the operators TL := —C + and Ti* := — C~ possess infinitely many 
zero-energy states, see Sect. 4.5. These are given by holomorphic and antiholomorphic 
functions for TC and TC*, respectively. In the absence of the coupling Hi = J2i A(gi, gi+i), 
a superparticle prepared in such a zero-energy state evolves in time without attenuation. 
Note that the unattenuated propagation of superparticles in holomorphic/antiholomorphic 
states in the bulk is strongly reminiscent of the semiclassical picture of electron guiding 
center drift along percolating equipotential lines. Our interest is in the critical long-range 
correlations. These will not be affected by the high-energy modes of TC and TC* - they 
do no more than give rise to complicated and irrelevant short-distance physics - but are 
determined by the zero-energy modes coupled together by the interaction Hamiltonian H\. 
Degenerate perturbation theory tells us how to construct the critical low-energy theory: we 
are to project the full space of quantum states of the a model on the alternating product 
of zero-energy modules V+i ® V-i <S> V+i <S> V-i <S> ■■■ • Calculating the matrix elements of Hi 
between the states of this multi-"spin" space, we will get the Hamiltonian of a quantum 
spin chain. 

To provide further motivation and background to Sect. 5.1, let us anticipate how 
the spin chain Hamiltonian is related to established theory. Starting from a quantum 
Hall system in the semiclassical high- field limit, where electrons drift slowly along the 
equipotential lines of a smooth random potential while performing rapid cyclotron motion, 
Chalker and Coddington || formulated a network model that mimics the relevant effects 
of quantum tunnelling near saddle points. The links of the network are directed, see Figure 
7a. Disorder is modelled by random 2x2 scattering matrices S connecting the incoming and 
outgoing channels at each node of the network. To bring the Chalker-Coddington model 
closer to the model that will emerge in Sect. 5.1, we take its Hamiltonian (or anisotropic) 
limit. This is done by setting S = exp(-ieh), with h some random hermitean 2x2 
matrix, and sending e —>■ 0. The resulting Hamiltonian dynamics is illustrated graphically 
in Figure 7b. Electrons run up or down on straight lines of alternating direction, and 
with a small probability amplitude (~ e) they tunnel from one line to a neighboring one, 
thereby reversing their direction of motion and acquiring a random change of phase. One 
can carry out the disorder average for this Hamiltonian version of the Chalker-Coddington 
model by using the supersymmetry method. Doing so, one arrives rather directly at the 
spin chain Hamiltonian of Sect. 5.1 [56], with electrons running up and down straight lines 
translating into degrees of freedom whose quantum states are specified by holomorphic 
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and antiholomorphic functions, respectively. 

The nonlinear a model is derived from the quantum limit of a random potential with 
correlation length l c short compared to the magnetic length Ib- We are discovering that 
its low-energy limit converges to a model associated with the opposite limit I c ^>Ib- This 
coincidence is a beautiful vindication of the hypothesis of low-energy universality at a 
critical point. 

5.1 Mapping on a Spin Chain 

We start with some basic definitions underlying the notion of what we are going to call a 
generalized "spin" chain. (We really ought to say "superspin" instead of just "spin" , but 
we prefer to use the shorter and less pretentious word when there is no risk of confusion.) 
As always, let E^ = 1, ...,p + q) be a set of canonical generators of End(C p+9 ) obeying 
the multiplication law EijE k i = 8j k Eu. End(C p+<3 ) is Z 2 -graded by \Eij\ = \i\ + \j\ mod 2, 
where \i\ = for 1 < i < p and \i\ — 1 for p + l<i<p + q. The introduction of a bracket 
operation ( "supercommutator" ) 

[Eij, E kl ] = 8 jk E u - (-)^ + ^ k ™5 H E kj (72) 

turns End(C p+9 ) into the Lie superalgebra gl(p, q; C), with the quadratic Laplace-Casimir 
element being Ylij(~ ) EijEji- gl(2n, 2n; C) can be identified with the complex Lie super- 
algebra of the symmetry group G of the quantum Hall model space G/K by the canonical 
isomorphism 

gl(2n, 2n; C) ^ gl(l, 1; C) ® gl(2) ® gl(n). (73) 

(Recall that the three factors in the tensor product relate to superspace, advanced-retarded 
space and extra space, in this order.) 

With V := V+i and the isomorphism (ff3[) being understood, we introduce E^ : V — > V 
by Eij — dT(Eij), where dT is the differential taken at the unit element of the represen- 
tation g i— > T g defined by Eq. fl5ip . The E^ are linear differential operators on the space 
of functions / G V. We shall refer to them as "spin operators" . They obey the super com- 
mutation relations 

[Eij y E k i\ = (-)M +mw+m 5j k Eu - 8 u E kj . (74) 

It is seen that these do not properly represent the superalgebra ([72]) , because of the appear- 
ance of overall extra minus signs. (Their origin is explained in Appendix C.) Nevertheless, 
the bracket (]7l| ) does define a Lie superalgebra and, as a matter of fact, it is a perfectly 
sensible object to work with. Operators E*j on the dual space V* := V_i are defined in iden- 
tical fashion. They, too, obey the superalgebra ([741). By Eq. ([46]), the set {£ij}i,j=i,...,4n 
naturally decomposes into four sets of operators A k i, B k i, C k i, and V k i (k,l = 1, ...,2n). 
The same decomposition is made for the spin operators E*j. For future reference we note 
the action of the spin operators on the special vectors v = V\ G V and v* = V\ G V*: 

djv = B* jV * = 0, 

A ij v = V ij v = A*jV*=V* j v* = (i^j), 

... 75 
A ll v = (-)^v/2, ViiV = 0, 

A*y = (-)i*i+V/2, v;y = o. 
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These express the fact that v and v* are lowest- and highest- weight vectors carrying the one- 
dimensional K-representations \i = exp STr n + In and = exp — STr tc + In, respectively. 

We turn to the formulation of the two-spin Hamiltonian. For two sites, the a model 
Hamiltonian H, Eq. (|3^), has V (g> V* for its degenerate space of strong-coupling ground 
states. The degeneracy is split by the two-site interaction Hi = a xx A(g, h)/a\ orthogonally 



projected onto V ®V* . Orthogonality is defined w.r.t. to the invariant scalar product (38), 
extended to tensor product space by 

(a®b*,c®d*) = (-)l b *H c l(a,c)(&*,cf). (76) 

We wish to express the projected Hamiltonian if e ff '■= PH\P in terms of spin operators. 
Since A(g, h) is G-invariant and invariance is preserved by orthogonal projection, H e n 
must be some Laplace-Casimir element of G formed from the spin operators on V and 
V*. The question is now, which Laplace-Casimir? Introducing Q(g) = gAg -1 we can 
write A(g, h) = STr Q(g)Q(h) /4. A symmetry transformation g \— > xg (x G G) takes Q(g) 
into xQ(g)x . The unique set of operators V — > V (V* — > V*) with this transformation 
behavior, which is identical to that of the adjoint representation, are the spin operators Sij 
It follows that if e ff must be proportional to the quadratic Laplace-Casimir element: 

H eS = jY t (-)W +1 £ ij ®s; i . 

ij 

This operator commutes with the two-spin generators Ski <8> 1 + 1 <8> as is easily verified 
by using the supercommutation relations ([74]) . By construction, i7 e ff is hermitean w.r.t. 
to the invariant scalar product (fFBl). (The interaction Hamiltonian A(g,h) obviously is, 
and hermitecity is preserved by orthogonal projection.) The constant of proportionality 
J = 4a xx /ai > is finite and is found by calculating some special nonzero matrix element 
of Hi. This is done in Appendix D. From now on we will often omit the tensor- pro duct 
symbol in multi-spin operators, keeping it for clarity however when writing multi-spin 
states. 

It is clear how to generalize all this to the case of 2N spins. The degenerate space 
of strong-coupling ground states becomes V ® V* <E> V <S> V* <g> ... (2A" factors), and the 
scalar product ( [76]) is extended in the natural way. (We just keep track of the minus signs 
that arise from changing the ordering of the fermions.) Orthogonal projection yields the 
multi-spin Hamiltonian 

Hs=Yl E(-) N+1 ( J +^i+i,jr + J-^Uji) , J+ = J- = J- (77) 

even I ij 

We have thus arrived at a spin chain reformulation of the nonlinear a 

model with topological coupling a* y = 1/2 (mod 1). The spin Hamiltonian (|77| ) is 
completely specified by the commutation relations ( [f4|) for the spin operators Sij and B*a 
and their action (|75| ) on the lowest- and highest- weight vectors v and v* . The spatial 
boundary conditions are chosen to be periodic. (In the case of open boundary conditions, 
we need to modify the modules V, V* for the edge spins in general; see Sect. 4.6.) According 
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to what was said at the end of Sect. 3.5, we can incorporate deviations of a xy from a* y by 
staggering the spin-spin interaction, i.e. by putting J± = (1 ± e) J with e ~ a xy — a* y . 

Because the spin chain is built from alternating modules V and V*, translation by 
one site cannot ever be a symmetry of the Hamiltonian Hs- However, for J_|_ = J_ and 
periodic boundary conditions, Hs possesses another symmetry which is just as good. There 
exists a canonical isomorphism T : V — *■ V*, f i— > / (complex conjugation). T squares to 
superparity by our using an adjoint of the second kind, / = (— y^f fl43[l . It extends to a 
map T : V ® V* <8> ••• — > V* <g>V ... in the natural way. From Eqs. (Q) and ( fi"5D one can 
show that the two-spin Hamiltonian behaves under complex conjugation as follows: 




This equation implies that the operator composed of translation by one site followed by 
complex conjugation, commutes with the Hamiltonian Hs of a periodic chain with J + = 
J_. 

The Hamiltonian Hs has been defined over the spin spaces V = V+\ and V* = V—\. 
Actually, we can formulate a whole family of Hamiltonians H s by taking V = V m and 
V* = V- m with m any positive integer. Since we know m to be the number of "edge 
states" , it is clear which role to assign to H s m ^ : it applies to a situation where m Landau 
levels are strongly mixed by the disorder. This is of experimental relevance for m = 2 since, 
for a weak magnetic field, the broadening caused by the disorder may exceed the Zeeman 
splitting and render Landau levels spin-degenerate. A plausible global phase diagram for 
this case has been suggested by D.K.K. Lee and Chalker |57j. (See also a recent preprint 



by D.H. Lee [|58| in this context.) Here we simply mention that, by using the coherent 
states of Sect. 4.5 to derive a functional integral representation of the m = 2 chain, we 
retrieve the nonlinear a model (13) with a xy = (mod 1). This implies a finite localization 



length in the center of the Landau band, in accordance with the conclusions reached in 



Ref. 57 



We have argued that projection of the a model on the superspin chain should retain the 
essential long wave length quantum Hall physics. Therefore, we expect the chain for m = 1 
and J+ = J_ where a derealization transition is known to occur, to be massless. Clearly, 
a direct proof of masslessness would be very welcome. A useful pool of ideas we can draw 



from is provided by the Lieb-Schultz-Mattis theorem [|59| for isotropic antiferromagnetic 
(ordinary) spin S = 1/2 chains, which has been generalized to arbitrary half-integer S by 
Affleck and Lieb |6(| . What these authors do is to construct a variational state (orthogonal 
to the ground state) with excitation energy of order 1/L for a chain of length L. By the 
variational principle, the existence of such a state implies that the Hamiltonian either has 
zero gap, or else a doubly degenerate ground state with a finite gap, in the thermodynamic 
limit L — > oo. Because of the resemblance between the S = 1/2 spin chain and the 
superspin chain with Hamiltonian (|77|), it seems that it should be possible to transcribe 
the theorem and its proof. Unfortunately, the transcription has not been done successfully 
up to now, for but one reason which however is vital: in order for the variational principle 
to apply, we need the Hamiltonian to be hermitean with respect to a quadratic form 
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which is positive definite. As we have seen, positivity is not a property enjoyed by the 
invariant scalar product (|38|). Nevertheless, as we shall see, the two-superspin Hamiltonian 
is completely diagonalizable and has positive excitation energies. We are thus encouraged 
to think that the difficulty is not of a fundamental kind but is a technical complication 
that can be overcome. 

To compute from the (super-)spin chain such physical observables as the conductance 
coefficients (g pq ), we need the spin chain expressions for the "contact operators" Oj of 
Sect. 2.2. For an ideal contact connected to a single spin, we know these to be given by 
O p = v(BnV, ■) for current, O q = Bi\v(v, ■) for voltage, and P v = v(v, ■) or P* = v*(v*, •) 
for a contact which is neither voltage nor current; see Sect. 4.6. In that subsection, we 
calculated the Hall conductance for a disk geometry using the independent edge spin 
approximation. Having constructed the Hamiltonian H$ governing the evolution of the 
interacting spin chain, we can now study the conductance coefficients in principle in the 
following, more interesting setting. We assume periodic boundary conditions in both time 
and space direction with lengths L and Lj, respectively. (Periodicity in space eliminates 
the boundary current.) Lq is chosen to be very large. We select two times ti < tf where all 
spins are connected to "particle reservoirs" . Recall that the effect of the reservoirs is to hit 
every degree of freedom in V or V* with a projector P v or P*, at the corresponding time. 
For Lq — ► oo, the time evolution outside the fixed interval [U,tf] projects on the ground 
state of H$- In Sect. 5.2 we will show that this ground state has zero energy and contains 
the state V'o = v ® v * ® v ® v * ® ••• with unit amplitude. To measure a conductance (g), 
we connect the k th F-spin, say, to a voltage probe at U and the I th V-spin to a current 
probe at tf. From the above expressions for O p and O q , we then get 

(g) = ^ f ,exp(-\t f -ti\H s )ipi) (78) 

where ipi (i/>/) is the multi-spin state obtained from ipo by creating one quantum of excita- 
tion v — > Bnv in the k th (I th ) l^-spin. When L\ is held fixed, (g) must decay exponentially 
with distance \tf — 1»| over some length £, the localization length. For the critical chain 
£ will be proportional to L\ for L\ large, (g) is expected to be governed by the laws of 
conformal invariance. However, since the definition of (g) involves two strings of operator 
insertions at the times t{ and tf, its behavior under conformal transformations is not ele- 
mentary. (The two strings impose boundary conditions.) Things simplify if we close off all 
the contacts which are neither voltage nor current. Then (g) turns into a spin-spin corre- 
lation function which measures conductance in a two-terminal (small contact) geometry of 
the kind discussed in Sect. 2.2. In this case, the hypothesis of conformal invariance relates 
the constant of proportionality between L\ and the localization length £ to the critical 
index governing the algebraic decay of (g) in the infinite system |§. 

To conclude this subsection, we mention that the a model expressions for the scal- 
ing fields of Sect. 2.1 are polynomials in the elements of the matrix- valued function 
Qij(g) = (gAg -1 )^; see e.g. Ref. [P5[ . Calculating the matrix element of such a polyno- 
mial Qi 1 j 1 Qi 2 j 2 -~Qirjr between two states in V or V*, we get an integral which is divergent 
for r > 2 because of the noncompactness of G/K. (The matrix elements between excited 
states of the monopole Hamiltonian are divergent for all r > 1.) In other words, these 
a model operators do not have a sensible restriction to the spin chain. (They are well- 
defined only when the system is perturbed by a finite frequency u -Me.) This is a technical 
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manifestation of the distinction between the different observables discussed in more direct 
physical terms in Sects. 2.1-2. 

5.2 Solution of the two-spin problem 

In Sect. 5.1 we constructed a spin Hamiltonian Hs by projection of the quantum Hamilto- 
nian of the nonlinear a model. Solving for the ground state and the low-energy excitations 
of Hs we will get an analytic description of the critical properties of integer quantum Hall 
systems. However, such a solution is not easy to find and has not been found yet. Here we 
will have to be satisfied with doing some preparatory work toward that ambitious goal. 

The two-spin Hamiltonian H := 1S an invariant second-order differ- 

ential operator acting on wave functions 

(/ <g> f*)(Z, Z, Z', Z') = v(Z, Z)v(Z', Z')<p(Z)ip*(Z') 

with holomorphic (p and antiholomorphic (p*. The following trick allows us to reduce the 
problem of diagonalizing if to a problem that has been solved before. Because / and /* 
are completely specified by the dependence of <p and tp* on Z and Z ', respectively, we 
may identify Z = Z' and Z = Z' without incurring any loss of information. On making 
this identification, H turns into an invariant second-order differential operator acting on 
functions 

F(Z,Z) = v(Z,Z) 2 <p(Z)<p*(Z). (79) 

There exists only one such invariant differential operator. This is the Laplacian C for 
G/K or, to be precise, C restricted to the truncated space of functions on G/K of the 
form ([F9|). Therefore, H = cC with some constant c and in the restricted sense just 
specified. In Appendix E we show c = —1/2. 

The Laplacian £ for Riemannian supersymmetric spaces G/K has been the object 
of some amount of study. A general procedure for the construction of its eigenfunctions, 
using such concepts as the Iwasawa decomposition of a semisimple complex Lie group 
and Harish Chandra's formula, has been described in ||48|1 - We are not going to review 
this procedure here but will specialize to n = 1 for brevity. In this case, an accidental 
simplification occurs which enables us to short cut the general theory. 

The key to making progress - for n = 1 as in general - is to consider the action 
of L on K-radial functions F(7i(g)) = F(-K{kg)) {k G K). By the Cartan decomposition 
] such functions can be regarded as functions on the positive part A + of 
a maximal abelian subgroup A C G with Lie algebra contained in V, the tangent space 
of G/K at 7r(l). In the present case, A ~ ]R xS 1 , and A + is conveniently parametrized 
by two coordinates x, y in the range 1 < x < oo and — 1 < y < 1. (These are the 
polar coordinates of Efetov, denoted by Ai, A in [^2|.) The origin 7r(l) has the coordinates 
x = y = 1. The invariant Berezin measure dgx on G/K induces a scalar product on the 
space of K-radial functions defined by 

dxjdy (x-y)- 2 F 1 (x,y)F 2 (x,y). (80) 



G = KA+K [50 
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With S being the function S(x,y) = x — y, the K-radial part Cft of C is the operator 

Cft = S (^£,^2 -\- Ci^j S 

where = 4d x (x 2 — l)d x and C^ 2 = 4<9 y (l — y 2 )d y are the radial parts of the Laplacians 
on the two-hyperboloid and the two-sphere, respectively. Note that Cft is hermitean w.r.t. 
the scalar product (|80| ) on the space of K-radial functions that vanish as x — y at the point 
x = y = 1. 

The separable form of Cft makes it possible to write down its eigenfunctions </>a, i 
immediately: 

<P\A x ^y) = ( x - y)^-(i\+i)/2(x)P(i-i)/ 2 (y)- 



Here, V v [x) is a Legendre function and P n (y) a Legendre polynomial |6l| . The eigenvalue 
is —(A 2 + I 2 ), and the quantum numbers have the range A G M + and / G 2^—1. Because 



of the special form of the functions ([F9D, not all of these quantum numbers occur in the 
spectral expansion of the two-spin Hamiltonian H. Since V = V+i and V* = V-i both 
carry angular momentum (or SU(2)-spin) < 1/2, the tensor product V <8> V* carries at 
most angular momentum 1. The largest angular momentum present in <p\,i is seen to be 
(I + l)/2. Therefore, the only /-quantum number that occurs for m = 1 is I = 1. (There is 
no restriction on A G 

The functions cf)\^ all vanish at the point x = y = 1 corresponding to Z = Z = 0. The 
two-spin quantum space V <S> V* contains the state v <S> v* whose wave function does not 
vanish at Z = Z = so, clearly, something is still missing. What is missing and must be 
added for completeness is the constant function <po = 1, which trivially is an eigenfunction 
of 

£# 4( ) 2 ( ^ ^ ~ ^ 9 ^ d 1 — y 2 d \ 
\dx (x — y) 2 dx dy (x — y) 2 dy ) 

with eigenvalue zero. Using <fio = 1 together with the functions <p\ i we can expand any 



K-radial function on G/K with "good" behavior at infinity. (See |I7| for more insight 
into the precise role played by <fio in a prototypical example.) 

Given a K-invariant stationary two-spin state ijj we can generate K-noninvariant sta- 
tionary states by applying symmetry transformations to the equation Htf; = Eip. Doing 
this in a systematic fashion by using the methods of supersymmetric Fourier analysis 
|£j],(|8| , we obtain a continuous set of G-multiplets labelled by A, / (with I = 1 being the 



only quantum number occurring for m = 1) that are eigenspaces of H and completely 
exhaust V (g> V*. 

To put this in perspective, recall the indefiniteness of the invariant scalar product 
(p8|), which jeopardizes the diagonalizability of H : V <g> V* — > V <8> V*. (In fact, it is 
precisely because of this indefiniteness that Laplace-Casimir elements of Lie superalgebras 
generically cannot be brought to diagonal form.) What saves the day is the existence of 
an auxiliary scalar product (•, on a restricted space (V <S> V*)# of K-invariant states. 
This auxiliary scalar product has all the attributes of a (Hilbert space) scalar product, and 
H : (V ® V*)^ — > (V <S> V*)^ is hermitean with respect to it. This in conjunction with 
G-symmetry guarantees the diagonalizability of H. 
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Note the pronounced "antiferromagnetic" character of H: it has a (G-) singlet ground 
state separated by a gap from a continuum of non-singlet excited states. The existence of 
a gap implies exponential decay of correlations in the two-spin system, which is related to 
exponential localization of all states in a Chalker-Coddington network model with only two 
links in the short, spatial direction. For the spin chain, exponential decay comes about 
because the repeated action of H on the state tpi = (B\\v) <S> v* creates on average an 
ever increasing number of quanta of excitation, eventually producing a time-evolved state 
which has exponentially small overlap with iftf = ipi. In the network model, exponential 
localization sets in when the elements of the transfer matrix become exponentially large. 
Thus, growth in the number of B- (and C*-) quanta is the spin chain equivalent of growth 
in the transfer matrix amplitudes of the network model. Using the Fourier decomposition 
of the Dirac delta distribution S D for G/K fl^] , 

+ 2* £ A *ff™ *», (81) 

I62M-1'' (A + ' > 

and doing the same calculations as in |[48|| , one can derive an exact expression for the 
(two-spin) conductance (g) of Eq. (f?8[). (g) decays as \tf — ti\~ 3 / 2 exp (— J\tf — ti\/2) for 
large distance \tf —ti\. 

Let us now present an independent argument showing why the ground state of H has 
zero energy. This argument has the merit of permitting easy generalization to many spins. 
We write the two-spin Hamiltonian in the form 

H = >J ) ' • 1 {AijAji + «yCi< + c a B ji + 

ij 

and consider the state v ® v*. It is a K-singlet since Tj~v = (i(k)v and T^v* = fi(k)~ 1 v*. 
Acting with H on it we get 

H(v ® v*) = ^(-) |l|+1 (By«) ® (C>*) , (82) 

ij 

as follows from Eqs. ([75]). Let Lm (M > 0) denote the (M + 1) -dimensional space of states 
generated by at most M actions of H on v ® v*. Since H commutes with the generators 
of K, all elements in Lm are K-invariant. Taking M — > oo we obtain a space, L^, which 
is closed under the action of H. (We here pay no attention to mathematical subtleties 
that might arise from the infinite-dimensionality of L^.) This space is decomposed by 
= Lm=o ® L c (direct and orthogonal sum). We claim that HL^ C L c . To prove this 
statement it suffices to prove e L c for %j) e L c , since H(y®v*) G L c from Eq. (^2|). Let 
Af = (An — T>n) /2 be the operator counting the number, N, of quanta of excitation 
in the first factor of V <8> V* . We have J\f(v ® v*) = 0. The commutation relations 
[Af, Bij] = Bij, [U.Cij] = -Cij, [N,Aij] = [A/\X>ij] = imply that H changes N by at 
most one unit. L c is a direct sum of spaces with N > 1. Let ijj now be a state with N = 1. 
There exists only one such state which has the additional property of being K-invariant. 
This is the state on the right-hand side of Eq. (|8~2D . Using the relations ([74]) and (|75j) , we 
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easily see that EijH 1 ^ 1 ^^ annihilates ip = Ey(-) |i|+1 (%^) ® so #V has 

zero intersection with v ® u*. But the image of any state ^eL c with more than one 
quantum of excitation in V carries at least one such quantum, which proves HL C C L c . 
We conclude that the cokernel of H : — > is (at least) one-dimensional. Therefore, if 
H can be brought to diagonal form - and we know it can be from the previous discussion 
- it must have an eigenstate with zero energy. On physical grounds, this eigenstate is 
expected to be the ground state. (Because the partition sum equals unity, the existence 
of a negative-energy state would give rise to unphysical spin correlations that increase 
exponentially with distance.) 

The result -f/Xoo C L c has the following consequence which is worth mentioning. Lq 
and L c are orthogonal spaces w.r.t. the scalar product ( |76|) . Hence, by the hermitecity of 
H, 

(H N (v ® v*),H N (v <8> v*)) = (v ® v*, H 2N (v <g> v*)) = 0, 

so the space L c is null, i.e. all of its elements have vanishing norm. 

The above argument readily generalizes to an arbitrary even number of spins. Let 
L denote the linear space spanned by i/jq = v Cg) v* ® ... ® v ® v*. If -0 is some state 
that differs from i/jq on more than two neighboring sites, then clearly Hgip fl Lq = {0}. If 
tp e Hf- X L Q (M > 1) differs from tp on no more than two neighboring sites, we are back 
to the two-spin problem analyzed before and again we have Hstp C\ Lq = {0}. This shows 
Hg 1 Lq to have zero intersection with Lq. Hence, the rank of Hs is reduced by one and we 
expect a zero-energy (ground) state as before. 

A neat application of the results of this subsection is the following. When the pro- 
cess of quantum tunnelling of electrons across saddle points is switched off, the Chalker- 
Coddington model (or, rather, a slight generalization thereof ||62|1 ) reduces to a network 
model which undergoes a classical percolation transition. This phase transition is in a dif- 
ferent universality class and has correlation length exponent v = 4/3 |k|. With a complete 
solution of the two-spin problem in hand, we can now write down a spin chain Hamiltonian 
for the classical percolation transition, too. To do so, it is easiest to make a checkerboard 
decomposition of the Hamiltonian and pass to a vertex model by exponentiation. The 
vertex model is specified by its R- matrix R : V <S> V* — > V* <S> V (Figure 8). Because of 
the connection of the supersymmetric vertex model with an underlying disordered network 
model, it is clear that the .R-matrix at the critical point of the percolation transition must 
be 

Rab*,c*d = ~ (8ad8b*c* + Sab*S c *d) • (83) 

(In the network model, electrons turn either right or left at a node, with equal probabilities. 
Correspondingly, the spin degrees of freedom of the vertex model must be transferred by the 
i?-matrix to the right or left, with equal weights.) On general grounds, the i?-matrix ( |83| ) 
must meet the requirement of G-invariance, and it does. With P : a®b* i— > (— )' & H a l&* eg) a 
the graded permutation operator, R can be written in the form R = P '(1 + Ho) /2 where Hq 
projects on the unique G-singlet state contained in the tensor product V ®V* . (Of course, 
this singlet state is nothing but the ground state of the two-spin Hamiltonian discussed 
above.) Restricting the resolution of unity ( |8~T| ) to V <S> V*, we obtain 



R = PR, R=l--[ d^(A)n ; 

2 Jr+ 
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Here dfi(X) = 2 2 Atanh(7rA/2)(A 2 + 1) 2 dX, and 11^ is the spectral projector onto total 
SU( 1,1) -spin with quantum number A, normalized by the condition 



([ dfi(x)h(x)n x 




[ d^(X)h(X) 2 U x 



To return to a Hamiltonian formulation, we have to compute the logarithm of the transfer 
matrix. We put R e = exp(—eH cp ) and get the two-spin Hamiltonian 



in the limit e — > 0. The Hamiltonian of the percolating chain is a translationally invariant 
sum of two-spin Hamiltonians. To move the system off the critical point, we stagger the 
interaction strength. 

The constancy of the excitation energy h(X) for H cp reflects the additional symmetries 
that distinguish classical percolation from the generic case of quantum percolation, where 
h(X) is nonconstant. We note that i/ qp = ) can be written in the form 



In this paper the study of the supersymmetric version of Pruisken's nonlinear a model for 
the integer quantum Hall effect has been resumed. By generalizing the work of Shankar and 
Read, we derived a Id Hamiltonian lattice formulation of the a model which has the right 
(naive) continuum limit. The lattice Hamiltonian consists of a site-diagonal part and a site- 
site interaction, the former being a supersymmetric generalization of the Hamiltonian for 
a charged particle moving in the field of a magnetic monopole. The monopole charge m is 
determined by the value of the Hall conductivity a xy . In the center of a Hall plateau region, 
where a xy is an integer, m vanishes for all sites in the bulk of the system and equals ±cr X y 
at the edges. In this special case we used a renormalization group assisted argument to 
calculate the conductance coefficient g pertaining to neighboring probes in a Corbino disk 
geometry. As is expected from the microscopic picture of directed transmission through 
channels localized near the sample boundary, we obtain either g = or g = \(J xy \ depending 
on the sign of a xy . This result suggests an interpretation of the monopole charge as the 
number of "edge states". 

For half-integral values of the Hall conductivity a xy corresponding to the points of 
transition between neighboring plateaus, the monopole charge in the bulk of the system 
alternates on sites (m = ±1). Projection onto the degenerate strong-coupling ground states 
of the lattice Hamiltonian yields a superspin chain Hamiltonian effectively describing the 
critical point of the a model. The space of quantum states for even and odd sites of 
the chain is a space of holomorphic and antiholomorphic functions V and V* , respectively. 
This alternating structure corresponds to the alternating orientation of links in the Chalker- 
Coddington network model, reflecting the directedness of guiding center motion in a strong 
magnetic field and a smooth random potential. A direct derivation of the superspin chain 
from the network model is possible ||56|1 . While Pruisken's nonlinear a model is derived 





(84) 



(ID with h(x) = a 2 + 1. 



6 Summary and Outlook 
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from the "quantum" limit of a random potential with correlation length l c short compared 
to the magnetic length lg, the network model applies to the opposite semiclassical limit 
l c Ib- Since both models map on the very same superspin chain, we conclude that 
the ratio I c /Ib is an irrelevant parameter for the IQHE transition. By the hypothesis of 
universality at a critical point, such irrelevance is not unexpected. 

We recalled the fact that the dynamic structure factor and other frequency-dependent 
correlation functions as well, have critical amplitudes that diverge at zero frequency. 
Within the field-theoretic framework this divergence is related to the fact that such cor- 
relators are not expressible as superspin correlation functions: the matrix elements of the 
corresponding a model operators between states in the superspin quantum spaces V and 
V* just do not exist. A quantity that does restrict to the superspin chain in a natural 
manner is d.c. conductance. 

While the nature of the subject is such as to force extensive use of mathematical 
argument, we tried to keep the mathematics to the minimum needed for clarity and com- 
pleteness. For example, we did not mention the role of the topological density as pullback 
of the Kahler form of the Kahler supermanifold G/K. The theory of connexions on her- 
mitian line bundles associated with a principal fibre bundle (central to Sects. 3.4-5) and 
of holomorphic sections of such bundles (central to Sects. 4.2-5), and the theory of Fourier 
analysis on supersymmetric spaces G/K (central to Sects. 4 and 5.2) were only mentioned 
parenthetically. We made these omissions in an effort to keep the paper within reasonable 
size and accessible to readers without a background in mathematical physics. 

Although the main body of this paper was concerned with the a model for 2d electrons 
in a strong magnetic field, most of what we did can be transcribed with minor modifica- 
tions to the symplectic a model describing time- reversal invariant systems with spin-orbit 
scattering. The role played by the state vectors v and v* is taken by the zero mode of the 



Laplacian discovered in Ref. [|48| . This zero mode is one state in an infinite-dimensional 
space of strong-coupling ground states of the quantum Hamiltonian of the symplectic a 
model. The strong-coupling ground states transform according to a representation of the 
symmetry group which is neither lowest-weight nor highest-weight nor irreducible. It is not 
clear at present whether the resulting superspin chain can serve as the universal critical 
theory of the Anderson metal-insulator transition for symplectic symmetry in d = 2. (The 
transition caused by spin-orbit scattering differs from that for quantum Hall systems in 
that a continuous symmetry is spontaneously broken.) 

An issue we briefly addressed is how the classical limit Ib — » fits into the superspin 
chain picture. This limit is not immediately accessible from our approach since the deriva- 
tion of Pruisken's model, as it stands, applies only to the limit Ib ^> l c - Nevertheless, 
by using the relation to the Chalker-Coddington model we were able to write down a su- 
perspin chain Hamiltonian for the classical percolation transition, too. This Hamiltonian 
has a higher degree of symmetry than the generic critical Hamiltonian describing quantum 
percolation. In this way, the superspin chain provides a unified framework in which to 
discuss both quantum and classical percolation and the crossover from the latter to the 
former. 

Since the mapping on a superspin chain retains the essential long wave length quantum 
Hall physics, the chain for a* y = 1/2 (mod 1), where a derealization transition occurs, 
should be critical. More precisely, we expect a unique G-singlet ground state with local 
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"antiferromagnetic" order and massless (non-Goldstone) excitations. The staggering of 
the spin-spin interaction for a xy 7^ o* xy favors singlets on strong bonds and opens a mass 
gap. Gaplessness of the critical chain should be guaranteed by the same mechanism that 
underlies the Lieb-Schultz-Mattis theorem for isotropic antiferromagnetic spin chains with 
half-integer spin. In spite of the tantalizing analogy, I have not yet succeeded in transcrib- 
ing the theorem and its proof. The problem is not lack of translational invariance but the 
indefmiteness of the invariant scalar product (j38[) , which prevents the application of the 
variational principle in its usual form. 

The ground-state energy of the superspin Hamiltonian vanishes identically for a chain 
of any length and for arbitrary boundary conditions, as is guaranteed by unbroken K- 
supersymmetry. The Virasoro algebra underlying the massless chain for m = 1 must 
therefore have central charge c = 0. Note that this does not imply triviality of a theory 
which is nonunitary. (The Kac-Moody algebra of conserved currents will have a nontrivial 
central extension.) 

In addition to the Hamiltonian (|77|) for m = 1, we can formulate quantum superspin 
chain Hamiltonians for all values of the monopole charge m > 0. By the interpretation of 
m as the number of transmitting channels, we know that such a generalized Hamiltonian 
describes a situation where m Landau levels are strongly mixed by the disorder. The chain 
for monopole charge m resembles a Id isotropic antiferromagnet with spin S = m/2. Large 
values of m offer the possibility of making a 1/m expansion. A natural starting point for 
such an expansion is the Neel state v®v* <g>v<g>v* Cg) ... . Approximating the spin operators 
by their large-m limits one gets a quadratic Hamiltonian which can be diagonalized by 
a Boboliubov transformation in the usual manner. The low-energy excitations of this 
Hamiltonian, the superspin waves, are found to have a linear dispersion relation, as one 
would expect from the close analogy to antiferromagnetic spin chains. These low-energy 
modes are the Goldstone modes restoring the continuous G-symmetry broken by the Neel 
state. Unfortunately, this type of analysis is not too useful for understanding the critical 
properties, for the 1/m expansion is ignorant of the expected profound difference between 
massive superspin chains for even m and massless ones for odd m. 

What are the prospects for obtaining some sort of exact solution? Although not too 
much work in this direction has been done as yet, a few statements can be made. The 
superspin chain Hamiltonian fl77|) is unlikely to be integrable. (The corresponding ordinary 
spin chains for spin S > 1/2 are not integrable.) It is therefore very fortunate that the 
universality we expect for m = 1 gives us much freedom to modify the superspin chain: any 
Hamiltonian which has the "right" symmetries and is not too remote in parameter space, 
should give the same low-energy physics — the physics of the unique RG-fixed point that 
describes the IQHE transition. Making use of this freedom, one may be able to construct 
a one-parameter family of commuting transfer matrices (with the superspin Hamiltonian 
being the logarithmic derivative at a special point) by solving the quantum Yang-Baxter 
equation (QYBE). This, unfortunately, is not immediately possible but requires an ex- 
tension of existing theory, for the following reason. We are looking for a G-invariant 
solution of QYBE. An efficient method for constructing such solutions has been described 
by Drinfel'd ||B5| , generalizing earlier work of Kulish et al. [^U . The basis of the method 



is an expansion of the .R-matrix R{u) : W\ (g) W2 — > W\ (g) W2 around u = 00 in the form 
R(u) = 1 + u~ x r + ... where u~ 1 r is a solution of the classical Yang-Baxter equation. The 
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crucial difference for our superspin chains is that no such classical limit exists. Because of 
the alternating structure of the chain, the .R-matrix is a map R{u) : V ® V* — > V* <8) V 
(Figure 8), which implies R{u) 7^ 1 for all u. (One might think that the isomorphism 
T : V — > V* , / 1— > / could be used instead of the identity. However, R(oo) = T fails to 
meet the requirement of G-invariance.) 

In my opinion, the most promising line of further attack is to study superspin chains 
with 1/r 2 exchange of the Haldane-Shastry type ]67j and/or explore the possible relation 



to a gauged WZW model. (Gauging is necessary to combat the indefinite metric on the 
noncompact symmetry group G.) 
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Appendix A: Proof of Formula ( |3gj ) 

Eq. (|33|) expresses C m by left-invariant vector fields. To switch to right-invariant vector 
fields we proceed as follows. In Eq. (|21|) we introduced a second-order differential operator 
X-p(d/dx) := X(d/dx) = g l i d 2 / dx^ dx l by inverting the invariant quadratic form B 
on V . Analogous differential operators Xg and A*; are defined by inverting B on Q = 
Lie(G) and K = Lie(K), respectively. These operators are related by the identity X-p = 
Xg — Xjc, which we insert into Eq. (|33D . By definition of the projection n : G — ► G/K, 
we have f(7r(gexp Y^x l ei)) = f(%(g)) for ^x l ei G /C. Moreover, fi(k(gexp ^afe^)) = 
(x(k(g))(x(expJ2 xle i) f° r Yy^ e i £ ^> an d 



X/c(d/dx)n(k(exp^2x' l ei)) 1 = Xjc(d/dx) exp (— m STr r K + ^x l ei) 

' x=0 ' 



= 

x=0 



by a simple calculation. Therefore, we may drop A^- Doing so we arrive at a formula for 
C m identical to (|33|) except that A = X-p is now replaced by Xg. After this preparation 
we can move exp ^x %e i to the left of g by using the invariance of Xg under exp ^2, x%e i 
g~ 1 (exp^2x l ei)g. This proves Eq. (|34]). 

Appendix B: Derivation of Eqs. (gl|), ([53]) and (|58|) 

We shall derive the coordinate expression ([51]) of the representation g 1— > T g , Eq. (|44|). 
Recall the definition of s : G/K — > Go = SU(n, n\2n) by 

i( A B \\\-( 1 BD^X ({I -BD- l CA- 1 )- 1 ' 2 A 

s{n{ {c Dj'l-yCA- 1 1 )\ (l-CA^BD- 1 )- 1 / 2 ) ' 

Solving the equation g = s(ir(g))k(g) for k, we get 

u( A b\ ((l-BD^CA-^A \ ( . 

*Hc d)>-\ (l-CA- 1 ^!)- 1 ) 1 ^ J • 
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From 7T + = (l + A)/2 we have fi(k) = exp(mSTr7r + Ink) = (SDet k) m / 2 exp (f STr A In k). 
For the second factor we obtain from (B.l) the simple formula 



exp^STrAlnA^ *))) = SDet(A/£>) 



m/2 



To express the factor SDet k(g), we take the superdeterminant of both sides of the equation 
g = s(irc(g))k(g), which gives SDetfc(^) = SDet g, since s(n(g)) G Go- Hence, 

/ A R \ / A R \ 

^(c J)) = SDet U d) SDet (^) m/2 - (B.2) 
From this we infer fj,(k(s(Z, Z))) = 1. Moreover, Eq. (B.2) in combination with 



A B\ 1 _ f 1 -A~ X B\ ( (A — BD~ X C)~ X 

C D J ~~ [-D^C 1 )[ (D-CA- 1 ^)- 1 



CB.3) 



yields fx(k(x)~ x ) = fx(k(x~ 1 )) for any x G G. Applying this identity to i = with 
/i = s(Z, Z), and using Eqs. (pESP, ( fI8p and (B.2) to calculate fi(k(s(Z, Z)~ 1 g)), we arrive 
at 

(\ m/2 

which proves Eq. ([Sip. 

The second calculation we do in this appendix is to work out the similarity trans- 
formation T+ = v^-T g v m for v m (Z,Z) = SDet(l - ZZ) m / 2 (to > 0). The effect 
of this transformation is to make the multiplier in Eq. (|5l] ) pick up an extra factor 
v m (Z, Z)~ 1 v m (g~ 1 ■ Z,g~ x ■ Z). To calculate it, we use the matrix function Q = gAg~ x = 
sAs~ x . From Eq. (D.l) it is seen that 1 — ZZ = ((1 + Q) bb/2) 1 where the indices BB 

mean the upper left (or Boson-Boson) block. Writing g = ( d) : ^~ 1= \C d'j anc ^ 
doing matrix multiplications, we find 

1 - GT 1 • Z){g~ x ■ Z) = Q(l + g- 1 Qg) B B^ =(A- ~ZC)~\\ - ZZ)(A' + B ' Z)~\ 

Now, from Eq. (B.3) we read off A' = (A - BD~ X C)- X and B' = -A'BD~ X , so taking the 
superdeterminant we obtain 

v m {Z,Z)- x v m {g- x ■ Z,g- X • Z) = SDet ( ^-C^BD- 1 ) ^ 2 
1 ' ; yy ' y J V (A- ZC)(D- ZB) ) 

Multiplying this with the multiplier in Eq. ( |5l|) and using SDet g = SDet(AD) SDet(l — 
CA- x BD~ l ) we get the multiplier in Eq. 

Finally, we derive the coordinate expression fl53|) of the monopole Hamiltonian 7i s = 
— £ m , Eq. (|33|). The pair of matrices X, X appearing in ( |53D is defined by the equation 
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exp(^x l ei) = s(D X XA,A 1 XD). Within the accuracy required by setting X = X = 
after differentiation, and because of the invariance of the differential operator d X x-> 
Eq. (13), under X ^ D~ l XA, X h-> A~ 1 XD, we may replace \{d/dx) in Eq. © by 
d x,x- Again, let g = ( c D J. Using Eqs. (H) we get 

g ■ (D^XA) = {Z + X)(1 + ZX)~\ 
g ■ (A^XD) = (Z + X)(1 + ZX)~\ 

where Z = CA~ X and Z = BD~ X are the coordinates of the coset 7r(g). This explains the 
argument of the function / on the right-hand side of Eq. (|53|) . We turn to the multiplier 
(x(k(g)k(g exp^x^i) -1 ) in Eq. (|33"1). The first factor, fx(k(g)), has been given in Eq. (B.2). 
For the second factor, we use 

/ A A- BX\ m / 2 
/i{k(gs(X,X))) = (SDet £) m / 2 SDet ' " ! ^ 



D + CX 



as follows easily from (B.2). Making the substitutions IhD 1 XA, Ini 1 XD and 
putting the factors together, we get the multiplier in Eq. ([53]). 



Appendix C: Representation of Lie superalgebras on super- 
functions 

We denote the superparity on any Z2-graded vector space V by | • |, i.e. \X\ = if X 6 V 
is even and |X| = 1 if X is odd. A complex Lie superalgebra Q is a Z2-graded vector space 
over C with a bracket operation satisfying the axioms 

\[X,Y] \ = \X\ + \Y\ mod 2, 

[x,y] = (-)W^i+ 1 [y,x], (c.i) 

[X, [Y, Z}] = [[X, Y},Z} + (-)\ X W Y \ [Y, [X, Z\\ 



An example is gl(p, q; C) defined by the bracket ([72]). Let now some Grassmann algebra 
A be given. Following Berezin ||43|| . we define the Grassmann envelope G(A.) of Q as the 
even part of the tensor product A <g> Q, that is to say, as the linear space of objects aX 
with a G A, X G Q and \a\ = \X\. When Q = gl(p, q;C), an element of ^(A) is what is 
commonly called a supermatrix of dimension (p + q) x (p + q) . The Grassmann envelope 
(?(A) can be turned E3] into a Lie algebra with Grassmann structure in either one of two 



ways: 

[aX, bY] := ab[X, Y], (1 st kind) 

[aX, bY] := a6(-) |x||b| [A,F]. (2 nd kind) 

For the definition of the second kind we imagine X and b to anticommute when both X and 
b are odd. The usual multiplication rules for supermatrices used in the physics literature 
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imply aXbY = abXY for a, b G A and X, Y G gl(p, C)- Therefore, the natural convention 
is to take the Grassmann envelope gl(p, q\A) to be a Lie algebra with Grassmann structure 
of the first kind. This is the convention we adopt. 

Now let Gc be a Lie group which has gl(p, q\A) for its Lie algebra, and let g t— > T g be 
a representation of Gc on some space of superfunctions. (For example, for p = q = 2n, Gc 
might be the complexified symmetry group of the quantum Hall model space G/K, and 
T might be the representation defined by Eq. (f44|).) Taking the differential at the group 
unit (t G R) 

dT(aX)f := —T expta xf 

at t=o 

we get a representation of the Lie algebra: 

[dT(aX), dT(bY)} = dT([aX, bY}). (C.2) 

dT(X) for X G gl(p, q; C) is defined by linearity, dT(aX) = adT(X). (This is the precise 
meaning of Eq. (p5|).) Thus, dT assigns to X a differential operator dT(X) acting on 
superfunctions /. Such an operator differentiates with respect to the commuting and 
anticommuting coordinates. If X is odd, so is dT(X). Consistency now requires us to 
observe the graded commutation law dT(X)b = (— )' x " 6 'kLT(AT) where b G A is some 
parameter. Hence, from (C.2) we deduce 

[dT(X),dT(Y)] = (-)\ x ^dT([X,Y]). 

It is easy to see that this bracket obeys the axioms (C.l), i.e. it defines a Lie superalgebra. 
Setting £ij = dT(Eij) with E{j the canonical generators of gl(p, q;C), we obtain the 
supercommutation relations [£ij,£ki] = (— )^ + ^^^ k ^ + ^ 5jk£u — 5uEkj, which are different 
from those for the Eij. We thus need to distinguish between two Lie superalgebras. The 
first one is generated by the E^, the second one by the £y. Our discussion clarifies how 
the two are connected: constructing from the former and the latter a Lie algebra with 
Grassmann structure of the first and second kind, respectively, we get Lie algebras that 
match, i.e. the latter represents the former. 

The sceptical reader might be reluctant to accept the violation of the representation 
property at the Lie superalgebra level. He is invited to convince himself that the only con- 
sistent way to fix the disease is to modify the basic law of multiplication of supermatrices. 
Clearly, such a modification would be very undesirable. 

Appendix D: Calculation of a Reduced Matrix Element 

We calculate the value of the exchange coupling constant J of the spin chain for m > 0. 
(Recall m = 1 for a spin-split and m = 2 for a spin-degenerate Landau level.) The 
matrix- valued function Q(g) = gAg~ x is expressed by 

(l + ZZ -2Z \fil-ZZ)- 1 

22 -i-zz)\ (l-ZZ)- 1 

in the coordinates Z, Z. We denote the matrix element (2Z(1 — ZZ)~ 1 )u by Qi\. From 
Sect. 5.1 we have the relation J = c^a xx ja\ where cq is the constant of proportionality in 
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PQi\v = cqBuv. (P projects onto V, and v = v m .) Making use of the invariant quadratic 
form (•,•) we get cq = (Biiv,Qi\v)/(BnV,Bnv). From Eq. (59a) for T+ = v~ 1 T g v we 
have (Bnv)(Z, Z) = mZuv(Z, Z). The invariant Berezin measure on G/K is expressed 
in the coordinates Z, Z by the flat measure dZdZ [12], so 

(B n v,Ql{v) = j dZdZ mZ n {2Z(l- ZZ)- x ) lx v{Z,Zf . (D.2) 

To calculate this integral, we use the identity 

m(Z(l - ZZ)- l ) xl v(Z, Zf = jv{Z, Z - ^u) 2 [ =o 

which follows from v (Z, Z) 2 = SDet(l — ZZ) m . By partial integration we get 

(B n v, Q\\v) = 2 J dZdZ v(Z, Zf = 2(v, v) = 2. (D.3) 

On the other hand, Eq. ( |69"D for a m+ 2,-m = Buv gives (Bnv, B\\v) = m. Hence, cq = 2/m 
and J = Aa xx /m 2 ai. 

The above calculation does not stand up to close scrutiny, for it turns out that the 
invariant Berezin integral in the coordinates Z, Z is plagued by boundary ambiguities |)8| 



in general. This means that the integrals in (D.2) and (D.3) must be corrected by the 
addition of extra terms. (I have analyzed the situation for n = 1 in some detail. In 
this case I find the correction terms for (D.2) to be nonvanishing only for m = 1,2. For 
m = 2, the correction term is precisely cancelled by the boundary term resulting from 
partial integration, so the first equality sign in (D.3) is correct as it stands in this case.) 



One can avoid dealing with this issue by switching to Efetov's polar coordinates |22| where 
the problem of boundary ambiguities is well understood [52,53]. Doing so, I find the final 
result to be correct for all m > 1 (and all n). 

Appendix E: Relation between the two-spin Hamiltonian and 
the Laplacian for G/K 



From Eq. (|3J), the monopole Hamiltonian Tt s = —C m is seen to have the expression 
H s = ^2ij(—)^ +1 £ij£ji, with the spin operators £^ = dT(Eij) being defined by the 
differential of T 3 , Eq. ([44]). Since H s is zero on its ground-state module V = V m , we may 
write the two-spin Hamiltonian H : V (g) V* —>■ V ®V* in the form 

# = !EH |l|+1 (^- + 

ij 

Returning to the notation used in Eq. fl34| ) we get 

(H(f ® /* )) (*(g)Mh)) = ~\\g{d / dx)ix(k(g)k(e^ e * g)' 1 ) 
x iiiKh^K^^h)) (f ® /*) (7r(e Exle ^), 7i(e^ xl ^h)) 

x=0 
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Now, following Sect. 5.2, we use the trick of identifying g = h. The multipliers then cancel 
each other, and switching from right- to left-invariant vector fields we arrive at the formula 

(H(f ® /*)) (Tr(g)) = -\\g{d/dx) (f ® /*) Me^N)) 

x=0 

= -\\ v {d/dx) (f ® /*) (^exp ^Ve*)) . 

x=0 

Comparison with formula ( pT|) proves that if coincides with minus one half times the 
Laplacian £ on G/K, after restriction to the subspace of functions which have the quantum 
numbers that are generated by the tensor product V £g> V*. 

Figure Captions 

Fig. 1: Three-terminal device (schematic) with small contacts. The configuration space 
Ai of the 2d electron gas is finite or infinite. The contacts Ci connect M. with particle 
reservoirs. 

Fig. 2: Lifting of a curve t i— > ^(t) on a principal fibre bundle G — * G/K with gauge 
connection (g~ 1 dg)ic- 

Fig. 3: Modification of the topological density. Aii is the union of the shaded regions. 
The topological coupling constant is doubled on Aii and set to zero on the empty regions. 

Fig. 4: Weight diagram of the lowest- weight module Vi- Weights are marked by full black 
dots. Weights enclosed by an additional circle are doubly degenerate. The lowest weight 
is (2, -2). 

Fig. 5: Corbino disk geometry with four terminals. Fat black lines symbolize G/K super- 
particles propagating in imaginary time. The particles at the outer and inner edges and 
in the bulk see monopole charges +m, —m, and 0, respectively. The sign of the monopole 
charge at the edges is indicated by arrows. 

Fig. 6: Ideal lead connected to a disordered quantum Hall sample. There are no direct 
transitions between incoming and outgoing channels. 

Fig. 7: (a) The network model of Chalker and Coddington. (b) Anisotropic limit of the 
network model. The dotted lines indicate tunnelling occurring with a small probability 
amplitude. 

Fig. 8: i?-matrix of the supersymmetric vertex model at the classical percolation transition. 
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